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ABSTRACT 


The  scope  of  this  thesis  involves  various  techniques  to  control  the  torque  of  a  voltage 
source  inverter-fed  induction  motor.  The  three  methods  investigated  here  are  voltage  by 
frequency,  known  as  “scalar  control”,  field-oriented  control  with  and  without  space 
vector  modulation,  and  direct  torque  control.  Although  all  three  techniques  have  been 
proven  effective,  each  technique  possesses  its  own  advantages  and  disadvantages.  In 
today’s  engineering  world,  time  and  money  are  saved  by  utilizing  software  developed  to 
accurately  model  a  physical  system  and  to  compare  with  hardware  previously 
implemented  in  a  lab.  Comparisons  can  be  made  to  determine  how  each  of  these  torque 
control  methods  perform  under  static  and  dynamic  modes  of  operation.  The  United  States 
military  can  benefit  from  this  study  by  modifying  or  adding  onto  the  software  models 
developed  here  with  little  cost  and  use  these  models  to  do  trade  studies  without  the  cost 
of  lab  demonstrations. 
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EXECUTIVE  SUMMARY 


The  Navy  is  constantly  looking  for  ways  to  improve  the  performance  and 
efficiency  of  shipboard  motors.  Due  to  the  current  state  of  our  national  deficit  and  the 
Navy’s  ever-shrinking  budget,  development  of  cheaper  alternatives  to  existing  motors 
continues  to  be  a  primary  focus.  These  alternatives  should  not  sacrifice  the  standards  for 
performance  set  by  machines  already  in  place.  If  we  can  maintain  or  improve  upon  these 
standards  by  making  them  more  robust  or  improving  their  dynamic  response  while 
eliminating  unnecessary  sensors  or  other  field  monitoring  devices,  the  money  that  can  be 
saved  will  be  substantial. 

The  evolution  of  DC  drives  to  various  forms  of  AC  drives  has  been  driven  by  the 
constant  need  for  improvement  in  performance,  simplicity,  and  reliability.  Although 
many  of  these  needs  appear  to  encompass  mutually  exclusive  goals,  new  developments 
strive  for  a  better  mix  of  advantages  versus  disadvantages.  A  simplistic  viewpoint  of  the 
evolution  of  drive  control  techniques  is  demonstrated  in  Figure  1. 


DC  Drive 


Flux  Vector  Control 


Figure  1.  Evolution  of  drive  control  techniques,  from  [1]. 
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The  goal  of  this  thesis  is  to  develop  an  induetion  machine  model  in  Simulink 
which  is  capable  of  accurately  implementing  torque  control  schemes  and  to  validate  this 
model  with  a  hardware  prototype.  It  is  possible  to  model  virtually  every  aspect  of  an 
induction  machine,  but  a  model  serves  only  as  a  means  to  reduce  time  spent  in  the 
laboratory  as  well  as  cost  and  risk.  A  useful  model  predicts  laboratory  results  and  circuit 
dynamics  but  can  only  be  declared  accurate  after  experimental  validation. 

The  construction  of  a  three-phase  voltage  source,  inverter-fed  (VSI-fed)  induction 
machine  drive  is  the  focus  of  this  thesis.  This  physical  system  was  tested  using  standard 
engineering  principles.  The  physical  parameters  of  this  machine  were  estimated  to 
develop  an  accurate  model.  A  typical  VSI-fed  induction  machine  drive  is  depicted  in 
Figure  2. 


Three-phase  Three-phase 

full  converter  inverter 


Figure  2.  Three-phase  VSI-fed  induction  machine,  from  [2]. 


The  objectives  of  this  research  were  to: 

•  Investigate  various  torque  control  schemes  for  an  induction  machine, 

•  Explain  the  control  principles  for  each  control  scheme  researched, 

•  Implement  a  physical  system  using  the  scalar  control  technique, 

•  Verify  a  Simulink  model’s  scalar  control  operation  by  comparing 
measured  results  from  the  hardware  prototype  to  the  simulated  results. 
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•  Compare  the  performance  of  various  torque  control  schemes  during 
transient  and  steady  state  conditions, 

•  Compare  speed,  current,  flux  and  torque  responses  of  field  oriented 
control  to  direct  torque  control, 

•  Provide  recommendations  as  to  which  control  technique  to  choose  based 
on  where  they  will  be  applied. 

To  control  speed,  torque,  and  position,  various  AC  drives  have  been  developed 
over  the  years.  Depending  on  quality  of  control,  these  systems  are  classified  as  either  low 
or  high  performance.  They  can  also  be  classified  as  scalar  (V/f)  or  vector-controlled 
based  on  their  control  principles.  In  general,  the  speed  of  an  AC  motor  depends  on  the 
frequency  of  the  stator  voltage,  while  the  developed  torque  is  typically  related  to  the 
stator  current. 

In  V/f  controlled  systems,  only  the  frequency  and  amplitude  of  the  fundamental 
stator  voltage  or  current  are  adjusted,  hindering  high  performance  of  the  system  under 
transient  operating  conditions.  Low  performance  V/f  drives  are  normally  employed  in 
machinery  such  as  adjustable  speed  pumps,  compressors,  and  fans  where  higher  quality 
control  is  unnecessary  [3]. 

The  main  concept  driving  V/f  control  is  to  keep  the  ratio  of  the  stator  voltage  to 
frequency  constant  to  maintain  constant  maximum  available  torque.  This  is  the  most 
popular  method  of  scalar  control  where  the  magnitude  of  the  variables  (frequency, 
voltage,  or  current)  is  controlled.  The  voltage  applied  command  is  calculated  directly 
from  the  applied  frequency  maintaining  a  constant  air  gap  flux  in  the  machine.  In  steady- 
state  operation,  the  flux  is  related  to  the  ratio  Vg/fg  where  I4  is  the  amplitude  of  the  motor 
phase  voltage  and  fg  is  the  synchronous  electrical  frequency  applied  to  the  motor.  This 
characteristic  is  defined  by  the  base  point  of  the  motor.  Below  this  base  point,  the  motor 
operates  at  optimum  excitation  because  of  the  constant  Vg/fg  ratio.  Above  this  point,  the 
motor  operates  under-excited  because  of  the  DC  bus  voltage  limit  [4]. 

Although  open  loop  V/f  control  was  utilized  in  this  thesis,  closed  loop  speed 
control  improves  performance  significantly.  With  closed  loop  control,  a  proportional- 
integral  controller  is  employed  to  improve  speed  accuracy  by  regulating  slip  speed  of  the 
motor  to  maintain  the  motor  speed  at  its  set  value. 


In  the  1990s  the  field  of  eontrolled  eleetrieal  drives  experieneed  rapid  expansion 
due  to  the  advanees  of  semieonduetors  in  power  eleetronics  eulminating  in  better 
mieroelectronie  mieroproeessors  and  digital  signal  proeessors.  These  technologieal 
advaneements  enabled  development  of  highly  effeetive  AC  drive  eontrols  with  lower 
power  dissipation  and  more  aceurate  eontrol  structures.  Drive  controls  became  more 
accurate  by  managing  the  three-phase  currents  and  voltages  with  “vector  controls.”  The 
most  common  form  of  vector  control  is  field-oriented  control  (FOC).  Three  major  points 
establish  this  scheme; 

•  Machine  current  and  voltage  space  vectors, 

•  Transformation  of  a  three-phase  speed  and  time-dependent  system  into  a 
two  coordinate  time-invariant  system, 

•  Effective  pulsewidth-modulation  (PWM). 

As  a  result  of  these  factors,  FOC  acquires  all  advantages  of  DC  machine  control 
and  frees  itself  from  the  mechanical  commutation  drawbacks.  Furthermore,  FOC  leads  to 
high  dynamic  performance  in  terms  of  response  times  and  power  conversion  by 
delivering  accurate  steady-state  and  transient  control.  A  basic  FOC  scheme  is  depicted  in 
Figure  3  [5]. 


Figure  3.  Basic  FOC  scheme,  from  [5]. 
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FOC  consists  of  controlling  the  stator  currents  represented  by  a  vector.  It  is  based 
on  projections  that  transform  a  three-phase  speed  and  time-dependent  system  into  a  two 
coordinate  {d  and  q  eoordinates)  time-invariant  system.  These  projections  led  to  a 
structure  similar  to  that  of  DC  machine  control.  FOC  machines  require  two  eonstants  as 
input  references:  the  torque  component  (aligned  with  the  ^-axis)  and  the  flux  component 
(aligned  with  the  J-axis).  As  FOC  is  based  on  projections,  the  control  structure  enables 
accurate  eontrol  in  all  operating  eonditions  by  seamlessly  handling  instantaneous 
eleetrical  quantities.  Thus,  FOC  solves  elassic  scheme  problems  by  easily  reaehing  a 
constant  reference  for  both  the  torque  and  flux  eomponents  of  the  stator  eurrent  [5]. 

Direct  torque  control  (DTC)  combines  FOC  theory,  direct  self-control  theory,  and 
advaneed  digital  signal  processing  (DSP)  and  applieation  speeific  integrated  circuit 
(ASIC)  technology  to  achieve  a  “sensorless”  variable-frequency  drive.  On  the  basis  of  the 
errors  between  the  reference  and  the  estimated  values  of  torque  and  flux,  it  is  possible  to 
direetly  control  the  inverter  states  in  order  to  reduee  the  torque  and  flux  errors  within 
prefixed  band  limits.  Unlike  FOC,  DTC  does  not  require  eurrent  regulators,  coordinate 
transformations  or  a  PWM  modulator  (as  a  eonsequence  timers  are  not  required).  In  spite 
of  its  simplicity,  DTC  allows  for  good  torque  control  in  both  steady-state  and  transient 
operating  conditions  [6].  A  basic  DTC  scheme  is  shown  in  Figure  4. 


Figure  4. 


Basie  DTC  scheme,  from  [6]. 


As  portrayed  in  Figure  4,  the  error  between  the  estimated  torque  T  and  the 
reference  torque  T*  is  the  input  of  a  three-level  hysteresis  comparator  and  the  error 
between  the  estimated  stator  flux  magnitude  (Pg ,  and  the  reference  stator  flux  magnitude 
is  the  input  of  a  two-level  hysteresis  comparator.  The  selection  of  the  appropriate 
voltage  vector  is  based  on  the  switching  table  given  in  Table  1. 

Table  1.  Basic  switching  table  for  DTC,  from  [6]. 


Sector 

1 

2 

3 

mm 

5 

6 

Cj  =  -1 

m 

V, 

h 

K 

Cj  —  +1 

k 

Vs 

Ct  =  ~1 

V, 

h 

Cj,  =  0 

K 

V, 

V, 

h 

F 

Crp  =  +1 

Vs 

V, 

V, 

Appropriate  use  of  the  optimum  switching  look-up  table  requires  knowledge  of 
the  stator  flux  linkage  space  vector,  because  the  sector  the  vector  resides  in  must  be 
known  [7].  This  simple  approach  allows  a  quick  torque  response  to  be  achieved; 
however,  the  steady-state  performance  is  characterized  by  undesired  ripple  in  current, 
flux  and  torque.  This  is  largely  due  to  the  absence  of  information  about  torque  and  rotor 
speed  values  in  the  voltage  vector  selection  algorithm  [6]. 

To  fairly  compare  FOC  to  DTC,  their  steady-state  and  step-response  for  three  step 
events  (at  t=0.5,  2,  and  3  seconds)  were  examined.  In  addition,  their  current,  torque,  and 
flux  responses  were  observed  closely.  These  responses  are  displayed  in  Figures  5-10. 
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Figure  5.  Speed  response  of  DTC  and  FOC. 
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Figure  6.  Phase  B  current  of  DTC  and  FOC  filtered  at  300  Hz. 
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Figure  7.  Reference  torque  for  DTC  and  FOC. 
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Figure  8.  Electrical  torque  for  DTC  and  FOC. 
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Figure  9.  Stator  flux  for  DTC  and  FOC. 
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Electrical  flux  for  DTC  and  FOC. 


Figure  10. 
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The  results  found  during  this  research  are  consistent  with  the  results  shown  in 
Table  2. 


Table  2.  Comparison  of  control  schemes,  from  [1], 


Control  Type 

Torque 

Control 

flux 

Control 

Response 

Advantages 

Disadvantages 

DC  Driv« 

Direct 

Direct 

High 

High  accuracy 

Good  torque  response 

Simple 

Motor  maintenance 

Motor  cost 

Encoder  required  for  high  accuracy 

Scalar  Frequency 
Control 

None 

None 

Low 

No  encoder 

Simple 

Low  accuracy 

Poor  torque  response 

Flux  Vector  Control 

Indirect 

Direct 

High 

High  accuracy 

Good  torque  response 

Encoder  always  required 

Direct  Torque 

Control 

Direct 

Direct 

High 

No  encoder 

Moderate  accuracy 

Encoder  required  for  high  accuracy 

Excellent  torque  response 


FOC  and  DTC  possess  distinct  advantages  over  V/f  control  and  suggest  that  scalar 
control  is  not  a  preferred  control  scheme  for  some  applications.  Although  FOC 
demonstrated  better  accuracy  (mainly  because  of  the  ability  to  control  the  current  with 
current  sensors),  this  accuracy  comes  at  the  expense  of  more  complicated  and  expensive 
electronics.  DTC  allows  for  direct  control  of  both  the  torque  and  the  flux,  but  it  is  not  as 
accurate  and  experiences  high  ripple  content  in  both  the  flux  and  torque  components.  The 
steady-state  current  of  DTC  was  significantly  lower  (less  electrical  losses)  than  that  of 
FOC  and  should  be  employed  with  machines  that  operate  at  constant  speeds  for  a 
majority  of  the  time.  FOC  would  better  serve  machines  that  require  frequent  speed 
changes  because  it  experiences  lower  peak-to-peak  currents  during  transients. 

It  was  demonstrated  that  by  using  simulation  to  verify  hardware  prototypes,  time 
and  money  can  be  saved.  To  further  demonstrate  this,  the  following  additional  research  is 
required: 

•  The  construction  of  a  prototype  that  utilizes  FOC  and  DTC, 

•  Performance  improvements  to  the  existing  model,  such  as  the  addition  of 
current  sensors  to  the  DTC  scheme, 

•  Application  of  this  model  to  machines  of  various  sizes. 

With  a  reliable  model,  the  Navy  has  a  means  to  build  a  large  scale  prototype  for 
laboratory  purposes  and  future  research. 
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I.  INTRODUCTION 


A.  BACKGROUND 

The  Navy  is  constantly  looking  for  ways  to  improve  the  performance  and 
efficiency  of  shipboard  motors.  Due  to  the  current  state  of  our  national  deficit  and  the 
Navy’s  ever-shrinking  budget,  development  of  cheaper  alternatives  to  existing  motors 
continues  to  be  a  primary  focus.  These  alternatives  should  not  sacrifice  the  standards  for 
performance  set  by  machines  already  in  place.  If  we  can  maintain  or  improve  upon  these 
standards  by  making  them  more  robust  or  improving  their  dynamic  response  while 
eliminating  unnecessary  sensors  or  other  field  monitoring  devices,  the  money  that  could 
be  saved  would  be  substantial. 

The  evolution  of  DC  drives  to  various  forms  of  AC  drives  has  been  driven  by  the 
constant  need  for  improvement  in  performance,  simplicity,  and  reliability.  Although 
many  of  these  needs  appear  to  encompass  mutually  exclusive  goals,  new  developments 
strive  for  a  better  mix  of  advantages  versus  disadvantages. 

1,  DC  Drive 

Torque  is  directly  proportional  to  armature  current  in  the  DC  motor.  Direct  torque 
control  becomes  possible  with  the  use  of  an  inner-current  loop.  Direct  flux  control  is 
possible  through  the  constant  magnetic  field  orientation  achieved  through  commutator 
action.  These  two  primary  factors  ensure  responsive  control  in  the  DC  drive.  Another 
advantage  of  the  DC  drive  is  that  only  simple  electronics  are  required  to  achieve  direct 
torque  and  flux  control.  On  the  other  hand,  initial  and  maintenance  costs  of  DC  motors 
are  high,  and  high  performance  speed  accuracy  can  only  be  achieved  by  including  an 
encoder  for  feedback  [1]. 

2.  Voltage-by-Frequency  (V/f)  Control 

Voltage-by- frequency  control  offers  the  advantage  of  operation  without  an 
encoder,  but  it  possesses  a  major  drawback  in  that  neither  torque  nor  flux  can  be 
controlled  either  directly  or  indirectly.  Control  is  provided  by  a  frequency  and  voltage 
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reference  generator  with  constant  Volts  per  Hertz  output.  This  generator  drives  a 
pulsewidth-modulated  (PWM)  modulator.  Unfortunately,  this  simple  arrangement 
provides  limited  speed  accuracy  and  poor  torque  response.  Flux  and  torque  levels  are  not 
controlled  by  the  drive  but  are  dictated  by  the  response  of  the  motor  to  the  applied 
frequency  and  voltage  [1]. 

3,  Field-Oriented  Control 

FOC  reestablishes  direct  flux  control  because  field  orientation  is  controlled 
electronically.  The  angular  position  of  rotor  flux  is  calculated  and  controlled  directly  by 
the  drive,  based  on  the  comparison  of  the  known  stator  field  vector  to  feedback  of  rotor 
angular  position  and  speed.  The  motor’s  electrical  characteristics  are  mathematically 
modeled  with  microprocessor  techniques  to  enable  processing  of  the  data.  However, 
torque  control  is  indirect  because  it  is  positioned  prior  to  the  vector  control  process  in  the 
control  algorithm. 

The  biggest  disadvantage  of  FOC  is  that  to  achieve  high-performance  speed  and 
torque  accuracy,  the  inclusion  of  a  pulse  encoder  is  necessary.  The  inclusion  of  the  PWM 
modulator,  which  processes  the  voltage  and  frequency  reference  outputs  of  the  vector 
control  stage,  creates  a  signal  delay  between  the  input  references  and  the  produced  stator 
voltage  vector.  These  last  two  factors  limit  the  ultimate  ability  of  FOC  to  achieve  rapid 
flux  and  torque  control  [1]. 

4,  Direct  Torque  Control 

Just  like  FOC,  DTC  reestablishes  direct  flux  control  and  implements  direct  torque 
control  as  well.  Both  flux  and  torque  are  controlled  by  hysteresis  controllers.  The  PWM 
modulator  stage  is  replaced  by  optimal  switching  logic,  and  thus,  the  delays  associated 
with  this  stage  are  removed.  The  original  benefits  associated  with  the  DC  drive  of  direct 
torque  control,  direct  flux  control,  and  high  responsiveness  are  all  reestablished.  As  a 
result,  better  torque  response  is  available  with  DTC.  Most  attractively,  if  moderate  speed 
accuracy  is  acceptable,  the  need  for  a  pulse  encoder  is  eliminated.  However,  DTC  has  its 
own  disadvantages  that  are  discussed  in  detail  later.  A  simplistic  viewpoint  of  the 

evolution  of  drive  control  techniques  is  demonstrated  in  Figure  1  [1]. 
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DC  Drive 


Scalar  Frequency  Control 


Flux  Vector  Control 


Direct  Torque  Control 


Figure  1.  Evolution  of  drive  eontrol  teehniques,  from  [1], 


B,  OBJECTIVE 

As  previously  stated,  there  are  many  sehemes  eapable  of  eontrolling  the  torque  of 
an  induetion  maehine.  The  purpose  of  this  thesis  is  not  to  make  eonelusions  as  to  whieh 
teehnique  is  the  most  effeetive  but  to  demonstrate  how  these  teehniques  eompare  both  in 
statie  and  dynamie  performanee. 

The  goal  of  this  thesis  is  to  develop  an  induetion  maehine  in  Simulink  capable  of 
accurately  implementing  torque  control  schemes  and  to  validate  this  model  with  a 
hardware  prototype.  It  is  possible  to  model  virtually  every  aspect  of  an  induction 
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machine,  and  a  model  serves  as  a  means  to  reduce  time  spent  in  the  lab  as  well  as  cost 
and  risk.  A  useful  model  predicts  laboratory  results  and  circuit  dynamics  but  can  only  be 
declared  accurate  after  experimental  validation. 

C.  APPROACH 

One  of  the  first  deeisions  to  be  made  was  to  determine  which  type  of  motor  would 
be  best  for  demonstrating  these  various  teehniques  of  torque  control.  One  of  the  machine 
theory  eourses  offered  at  the  Naval  Postgraduate  School  examined  a  three-phase  Voltage 
Source  Inverter-fed  (VSI-fed)  induetion  machine  and  determined  its  cireuit  parameters.  A 
schematie  of  this  machine  is  shown  in  Figure  2  and  is  discussed  in  further  detail  in 
Chapter  II. 


Three-phase  Three-phase 

full  converter  inverter 


After  ealeulating  this  maehine’s  parameters,  it  was  modeled  in  Simulink.  Due  to 
the  scope  of  this  thesis,  only  V/f  control  was  implemented  in  the  laboratory.  Upon 
verifying  that  simulated  and  measured  results  utilizing  V/f  control  matched  elosely,  we 
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assumed  this  model  behaved  aecurately  under  all  torque  control  schemes.  This  model  was 
used  to  simulate  the  other  control  schemes  (i.e.,  FOC  and  DTC),  and  all  three  simulated 
torque  control  methods  were  compared. 

D,  THESIS  ORGANIZATION 

Induction  machine  theory  and  how  the  Simulink  model  was  developed  are 
presented  in  Chapter  II.  Theory  behind  the  various  torque  control  techniques  are 
addressed  in  Chapter  III.  The  operational  modes,  advantages  and  disadvantages  of  each 
are  explained  in  detail.  Proper  functionality  of  the  Simulink  model  under  V/f  control  is 
verified  through  the  comparison  of  simulated  and  experimental  data  and  discussed  in 
Chapter  IV.  The  V/f  simulation  is  used  to  demonstrate  that  both  the  FOC  and  DTC 
simulations  accurately  predict  the  behavior  of  this  induction  machine.  A  comprehensive 
comparison  between  FOC  and  DTC  is  included  in  Chapter  IV.  Conclusions  and 
recommendations  for  future  research  are  discussed  in  Chapter  V.  Data  sheets  for  the 
equipment  used  and  the  Matlab  code  needed  to  regenerate  the  results  found  here  are 
included  in  the  appendices. 
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II.  INDUCTION  MACHINE  THEORY 


A.  INTRODUCTION 

Simulink  is  a  powerful  tool  that  can  be  used  to  simulate  any  system  that  can  be 
represented  mathematically.  To  understand  the  simulated  torque  control  techniques 
presented  here,  it  is  necessary  to  first  possess  a  thorough  understanding  of  the  equations 
that  represent  the  induction  machine.  The  equations  used  to  characterize  this  machine  are 
discussed  in  further  detail  in  [2].  Presented  here  are  the  induction  machine  voltage 
equations,  flux  linkage  relationships,  as  well  as  the  equations  for  stator,  rotor,  and  mutual 
winding  inductances.  When  creating  an  equivalent  circuit,  a  standard  convention  is  to 
refer  the  rotor  variables  to  the  stator  windings  using  the  appropriate  turns  ratio.  This 
convention  applies  to  the  magnetizing  and  mutual  inductances  as  well.  More  can  be  read 
concerning  induction  machine  theory  as  used  in  this  thesis  application  in  [3]. 

1.  From  DC  to  AC 

In  the  past,  DC  motors  were  used  widely  in  variable-speed  applications  because 
their  flux  and  torque  can  be  controlled  by  the  field  and  armature  current.  The  existence  of 
the  commutator  and  brushes  presents  many  challenges  for  DC  machines  because  these 
parts  require  periodic  maintenance  and  can  only  be  used  in  certain  environments.  Their 
limitations  under  high  speed,  high  voltage  operational  conditions,  however,  can  be 
overcome  by  the  application  of  AC  motors.  AC  motors  are  robust,  unaffected  by  heavy 
overloading,  and  easy  to  maintain  while  being  relatively  inexpensive.  They  can  be 
designed  with  substantially  higher  output  ratings  for  low  weight  and  low  rotating  mass 
based  on  their  size  alone  [4]. 

2,  Early  Misconceptions 

Torque  in  AC  and  DC  machines  is  similarly  produced.  This  was  not  emphasized 
before  the  1970s.  As  a  result  FOC  did  not  materialize  earlier.  The  precept  given  in  many 
popular  textbooks  on  machine  theory  falsely  implies  that  the  monitoring  of  the 
instantaneous  electromagnetic  torque  of  an  induction  machine  required  monitoring  of  the 
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rotor  currents  and  position.  Even  more  reeent  publieations  substantiated  these 
misconceptions  primarily  because  the  convoluted  expression  for  the  instantaneous 
electromagnetic  torque  had  not  yet  been  simplified.  Fundamental  physical  laws  and 
space-vector  theory  proved  that  torque  eould  be  expressed  as  the  product  of  a  flux 
producing  current  and  a  torque  produeing  current  if  a  flux-oriented  referenee  frame  was 
utilized.  In  this  case,  the  stator  current  components  (normally  expressed  in  a  stationary 
reference  frame)  are  transformed  into  a  new  rotating  reference  frame  which  rotates  along 
with  a  seleeted  flux-linkage  spaee  vector  [4]. 

3.  FOC  and  DTC 

The  torque  in  both  DC  and  AC  maehines  is  controlled  by  the  motor  eurrents. 
Unlike  DC  maehines,  both  the  phase  angle  and  the  modulus  of  the  eurrent  must  be 
controlled  in  AC  machines.  This  brings  reason  to  the  term  “vector  control.”  In  addition, 
where  the  commutator  and  brushes  of  a  DC  maehine  fix  the  orientation  of  the  field  flux 
and  armature  magnetomotive  force  (analogous  to  emf),  in  AC  maehines,  the  field  flux 
and  spatial  angle  of  the  armature  magnetomotive  force  require  external  control.  Without 
this  external  eontrol,  the  spatial  angles  between  the  various  fields  in  an  AC  maehine  vary 
with  the  load  and  ereate  unwanted  oscillating  dynamic  response.  Through  external 
control,  the  torque  and  flux  producing  current  components  are  deeoupled,  and  the 
transient  response  characteristies  emulate  those  of  a  separately  excited  DC  maehine.  This 
system  will  adapt  to  load  disturbanees  and/or  reference  value  variations  just  as  quiekly  as 
a  DC  machine. 

Vector  control  is  only  one  type  of  torque  eontrol.  The  other  type  of  high 
performanee  torque  control  is  DTC.  DTC  is  achieved  by  direct  and  independent  control 
of  the  flux  linkages  and  electromagnetie  torque  through  the  selection  of  optimal  inverter 
switehing  modes,  whieh  gives  fast  torque  response,  low  inverter  switehing  frequency,  and 
low  harmonic  losses.  How  these  control  methods  compare  with  V/f  control  and  also  with 
one  another  in  both  transient  and  dynamie  operating  eonditions  is  demonstrated  in  this 
thesis  [4]. 
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B, 


VSI-FED  INDUCTION  MACHINE 


The  model  used  for  this  thesis  and  how  its  parameters  were  calculated  are 
described  in  this  section.  This  VSI-fed  induction  machine  was  implemented  in  the 
laboratory,  and  an  equivalent  model  was  designed  using  Simulink.  After  verifying  V/f 
control  of  this  machine  through  measured  and  simulated  results,  we  assumed  that  the 
model  behaves  accurately  regardless  of  which  torque  control  scheme  is  utilized.  Then,  the 
Simulink  model  was  modified  to  simulate  FOC  and  DTC  separately  and  then  compared 
with  V/f  control  for  accuracy. 

1.  Balanced  Three-Phase  Machines 

The  three-phase  induction  machine  has  three-phase  windings  distributed  on  the 
stator  iron  and  either  three-phase  windings  wound  on  the  rotor  or  an  equivalent  squirrel- 
cage  structure.  Applying  balanced  three-phase  power  to  a  balanced  AC  motor  implies 
that  the  line-to-line  voltages  across  the  motor  terminals  are  sinusoidal  in  steady-state.  In  a 
balanced  AC  motor,  it  can  be  shown  that  the  motor  voltages  are  also  sinusoidal  and, 
similarly  to  the  motor  currents,  the  three-phase  voltages  sum  to  zero  on  an  instantaneous 
basis  as  expressed  by 

^an(t)  +  +  ^^cn(t)  =  0  (1) 

where  and  are  these  motor  voltages. 

This  is  also  valid  on  an  average  basis 

\n{^)  +  Vl,„{t)  +  Vcn{^)  =  ^  (2) 

when 

4  (0  +  4  (0  +  4  (0  =  0  (3) 

where  4,4  4  are  the  phase  currents  [5]. 

The  balanced  power  applied  to  the  stator  circuits  produces  a  rotating  magnetic 
field  (North/South  pole  pair)  which  cuts  the  windings  on  the  rotor  and  induces  voltages. 
The  windings  are  typically  short  circuited  and  current  flows.  In  turn,  the  rotor  currents  set 
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up  a  rotating  magnetic  field  relative  to  the  rotor.  The  interaction  of  these  fields  is  the 
mechanism  of  torque  production  in  induction  machines  [6]. 

2,  The  J-Equivalent  Circuit 

One  of  the  machine  theory  courses  offered  at  Naval  Postgraduate  School  required 
students  to  perform  a  laboratory  experiment  designed  to  calculate  the  parameters  of  a 
VSI-fed  induction  machine.  From  a  simulation  and  analysis  perspective,  it  is  desirable  to 
calculate  the  approximate  parameters  for  the  induction  machine  T-equivalent  circuit.  To 
do  so,  no-load  and  locked-rotor  tests  were  performed.  The  T-equivalent  circuit  for  steady- 
state  operation  of  a  symmetrical  induction  machine  is  portrayed  in  Figure  3. 


Figure  3.  Equivalent  circuit  for  steady-state  operation  of  a  symmetrical  induction  machine, 

after  [2]. 


When  an  induction  machine  runs  at  no-load,  the  slip  is  exceedingly  small  which 
implies  that  the  equivalent  resistance  is  very  large.  As  a  result,  the  rotor  current  is 
negligible,  and  the  equivalent  circuit  reduces  to  only  the  magnetizing  branch  which  is 
modeled  by  a  reactance  and  a  parallel  resistance  R^,  in  series  with  the  stator 

resistance  and  leakage  reactance.  The  resistance  allows  the  core  losses  of  the  machine 

to  be  taken  into  account.  The  values  of  this  reduced  circuit  were  derived  by  applying 
rated  line-to-line  voltage  and  measuring  the  no-load  current  and  total  three-phase  active 
power. 
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When  the  rotor  of  an  induction  motor  is  locked,  the  stator  current  is  on  the  order 
of  six  times  greater  than  the  rated  current.  Furthermore,  the  slip  is  equal  to  one.  While  the 
slip  is  equal  to  one,  the  resistance  is  small  compared  to  the  magnetizing  branch.  Thus,  the 
magnetizing  branch  was  neglected. 

3,  Rated  Values  and  Average  Phase  Resistance 

The  rated  current  and  rated  phase  voltage  of  the  stator  winding  as  well  as 
the  machines  rated  speed,  electrical  frequency,  and  horsepower  of  the  motor  are  listed  in 
Table  1. 


Table  1.  Rated  values. 


Parameter  name 

Value 

Rated  current  (stator 
winding) 

1.5  A 

Rated  phase  voltage 
(stator  winding) 

120Vnm 

Rated  speed,  0)^ 

1670  rpm 

Rated  electrical 
frequency', 

60  Hz 

Horsepower 

14  HP 

First,  the  number  of  poles  was  calculated  from 
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p  =  {2fx60)^NLS  (4) 

where  p  is  the  number  of  poles,  is  the  eleetrieal  frequency  in  Hertz,  and  NLS  the  no- 

load  speed  was  measured  in  rpms.  With  a  no-load  speed  of  1789  rpm,  it  was  discovered 
that  this  was  a  four  pole  machine,  from  the  horsepower  and  speed,  the  torque  can  also  be 
calculated  from 

P  =  m^T  =  —  (5) 


where  P  is  %  horsepower,  T  is  the  calculated  torque  and  co^  is  the  rated  speed  in 
radians/second.  This  calculation  leads  to  a  torque  of  9.4  Ibf  •  in. 

By  measuring  the  resistance  of  each  stator  phase  the  average  stator 
winding  resistance  can  be  determined  from 

n  _  ^Ph\  +  ^ph2  +  Rphl, 

^ph  ~  2 

where  is  the  average  phase  resistance  in  ohms  and  is  1 1.91  Q. 


a.  No-Load  Measurements 

The  circuit  used  to  measure  no-load  characteristics  is  illustrated  in  Figure 


Figure  4.  Circuit  for  measuring  no-load  characteristics,  from  [6]. 
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By  using  the  expression  for  instantaneous  three-phase  power  and  the  fact 
that  the  machine  is  wye-connected,  the  setup  in  Figure  3  yields  total  steady-state  three- 
phase  power 


'■^ashs+^bshs+^cshs- 


(7) 


This  can  be  manipulated  to  obtain 


(8) 


where  ,  and  are  the  stator  phase  voltages  and  4*  >  hs  >  4*  are  the  stator 

phase  currents. 

By  applying  200  VAC  and  recording  the  line-to-line  voltages,  line 
currents,  power  readings,  and  machine  speed  the  total  three-phase  power  was  calculated 
to  be  138  W  (no-load  characteristic).  These  no-load  measurements  are  listed  in  Table  2. 


Table  2.  No-load  measurements. 


Parameter  name 

\’alue 

Channel  A  Line  Current 

0.684  A 

Channel  A  Power 

60  W 

Channel  B  Line -to -Line 
Voltage 

201  V 

Channel  B  Line  Current 

0  6S  A 

Channel  B  Power 

78  W 

Machine  Speed 

1789rpm 
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b.  Locked-Rotor  Measurements 

For  the  locked-rotor  test,  the  line-to-line  voltage  was  adjusted  to 
approximately  60  VAC,  and  120  VAC  was  supplied  to  the  dynamometer.  The  motor  was 
connected  to  the  dynamometer  with  a  timing  belt,  and  the  dynamometer  was  set  to 
provide  maximum  load.  After  measuring  the  appropriate  currents,  powers,  and  voltages, 
the  total  three-phase  power  was  computed  to  be  234  W  (locked-rotor  characteristic). 
These  locked-rotor  measurements  are  listed  in  Table  3. 

Table  3.  Locked-rotor  measurements 


Parameter  name 

Value 

Channel  A  Line-to-Line 
V’oltage 

99  V 

Channel  A  Line  Current 

2.09  A 

Channel  A  Power 

42  W 

Channel  C  Line-to-Line 

V’  oltage 

100  V 

Channel  C  Line  Ciurent 

2.04  A 

Channel  C  Power 

192  W 

c.  Final  T-Equivalent  Circuit  Parameters 

From  the  data  collected,  the  following  parameters  for  the  T-equivalent 
induction  machine  circuit  were  calculated;  rs.Xis.X'i^.X^  and  the  core  loss  resistance  . 
The  rms  phase  A  stator  voltage  v^  was  determined  from 


V 

V 

as  ^3 


(9) 
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where  is  the  rms  line  a-to-line  b  voltage  and  was  116  Vrms- 


To  calculate  the  no-load  power  of  the  core  PaoreNi^  the  following 
expression  was  used 


P  =P  -3  1  r 

^  core,NL  ^  in,NL  -’^as,NL  's 


(10) 


where  P.^  is  the  no-load  input  power,  is  the  no-load  phase  A  stator  current,  and 

is  the  resistance  of  the  stator.  A  no-load  power  of  the  core  of  121.5  W  was  computed. 
Next,  the  core  resistance  i?.,was  determined  from 


R  =■ 


3 


V, 


as.NL 


(11) 


core.NL 


where  is  the  instantaneous  no-load  phase  A  stator  voltage.  This  results  in  a  core 

resistance  of  332  Q. 

The  phase  angle  difference  AO  was  calculated  using  the  rated  frequency 

of  60  Hz  from 


360°At 


// 


(12) 


where  At  is  the  measured  difference  between  the  voltage  and  current  phase  (3.92  msec). 
This  phase  angle  difference  is  84.7°. 


Starting  with  the  expression 


=  ^r/+(X^  +  X,f 


(13) 


where  is  the  mutual  reactance  and  is  the  stator  leakage  reactance,  and 

rearranging  Equation  (13)  in  the  following  way 


Y 

^  as,NL 

^  as,  NL 

(14) 
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we  can  determine  the  total  reactance.  A  total  reactance  of  170  Q  was  calculated.  In  turn, 
the  resistance  of  the  rotor  on  the  secondary  side  can  be  found  using  the  expression 


s 


P:.. 


—  r. 


3  / 


,LR 


(15) 


where  is  the  locked-rotor  input  power  and  is  the  locked-rotor  instantaneous 

phase  A  stator  current.  This  computation  resulted  in  a  resistance  of  5.95  Q.  Rearranging 
Equation  (13)  into  the  form 


V 

^  as,LR 

- 

2 

^ as, LR 

M 

V  ^  J 

(16) 


and  Equation  (16)  into 


E 


as,LR 


s 


+ 


J 


(K+x„)" 


(17) 


where  X  is  the  rotor  leakage  reactance  on  the  primary  side  and  is  the  stator 
leakage  reactance  (they  are  equal)  makes  calculating  the  reactance  possible.  Rearranging 
Equation  (17)  allows  for  both  to  be  solved  from 


V 

^  as,LR 

p 

2 

f 

lE 

^as,LR 

"  1 

V 

J 

1 

7 

(18) 


and  X\^  and  X^^  are  equal  to  10  Q  apiece.  Einally,  the  mutual  reactance  Xj^  can  be 
obtained  as 


(19) 


Given  that  Xf^  is  10  Q  it  is  easy  to  see  that  there  is  a  total  mutual  reactance  of  160  Q. 
The  calculations  performed  above  are  summarized  in  Table  4. 
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Table  4.  Circuit  parameter  values. 


Parameter  name 

Value 

Rotor  resistance, 

5.95  Q 

Stator  resistance, 

5.95  Q 

Stator  linkage  reactance, 

10  Q 

Mutual  reactance, 

160  Q 

Rotor  linkage  reactance,  X,,. 

10  Q 

This  is  the  T-equivalent  circuit  with  coil  one  selected  as  the  reference  coil.  The  values  on 
the  rotor  side  ( and  )  have  been  scaled  by  the  turns  ratio  [6]. 

Once  the  values  of  the  T-equivalent  circuit  were  calculated,  these  values 
and  the  specifications  of  the  machine  used  in  this  lab  were  used  to  develop  an  accurate 
model  of  the  machine  in  Simulink.  To  ensure  that  the  model  accurately  predicted  the 
behavior  of  this  machine,  these  same  parameters  were  incorporated  into  the  model.  The 
simulated  and  measured  results  were  then  compared  to  ensure  that  the  model  accurately 
predicted  the  behavior  of  the  physical  machine.  These  results  are  shown  in  Figure  5. 


Figure  5.  Measured  and  equivalent  electrical  torque  versus  speed. 
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III.  TORQUE  CONTROL  METHODS  FOR  VSI-FED  INDUCTION 

MACHINE 


A,  VOLTAGE-BY-FREQUENCY  CONTROL 

To  control  speed,  torque,  and  position,  various  AC  drives  have  been  developed 
over  the  years.  Depending  on  quality  of  control  these  systems  are  elassified  as  low  or 
high  performanee.  They  ean  also  be  elassified  as  scalar  (V/f)  or  vector-controlled  (FOC) 
based  on  their  control  principles.  In  general,  the  speed  of  an  AC  motor  depends  on  the 
frequency  of  the  stator  voltage,  while  the  developed  torque  is  typieally  related  to  the 
stator  eurrent. 

In  V/f  eontrolled  systems,  only  the  frequeney  and  amplitude  of  the  fundamental 
stator  voltage  or  current  are  adjusted,  hindering  high  performanee  of  the  system  under 
transient  operating  eonditions.  These  low  performance  V/f  drives  are  normally  employed 
in  machinery  sueh  as  adjustable  speed  pumps,  compressors,  and  fans  where  higher 
quality  eontrol  is  unnecessary  [5]. 

1.  V/f  Concept 

The  ratio  of  the  stator  voltage  to  frequency  is  kept  constant  to  maintain  maximum 
available  torque  in  variable-speed  induetion  motor  drives.  This  is  known  as  V/f  eontrol 
and  is  the  most  popular  method  of  sealar  control  where  the  magnitude  of  the  variables 
(frequency,  voltage,  or  current)  is  controlled.  The  voltage  applied  command  is  calculated 
direetly  from  the  applied  frequeney  maintaining  a  constant  air  gap  flux  in  the  maehine.  In 
steady-state  operation,  the  flux  is  related  to  the  ratio  Vg/fg  where  I4  is  the  amplitude  of 
the  motor  phase  voltage  and  fg  is  the  synchronous  electrical  frequency  applied  to  the 
motor.  This  oharaeteristie  is  defined  by  the  base  point  of  the  motor.  Below  this  base 
point,  the  motor  operates  at  optimum  excitation  because  of  the  constant  Vg/fg  ratio. 
Above  this  point,  the  motor  operates  under  exeited  beeause  of  the  DC  bus  voltage  limit 

m. 
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a.  Faraday  Law  and  the  Speed  of  an  Induction  Motor 

The  speed  of  an  induetion  motor  can  be  controlled  by  varying  the 
frequency  of  the  three-phase  power  supply.  To  maintain  constant  flux  density,  the  applied 
voltage  must  be  changed  in  the  same  proportion  as  the  frequency  according  to  Faraday’s 
Law.  Any  change  in  the  magnetic  environment  of  a  coil  or  wire  will  cause  a  voltage 
(emf)  to  be  induced  in  the  coil.  Although  there  are  many  ways  to  change  this 
environment,  this  voltage  is  generated  regardless  of  what  causes  this  change.  The 
magnetic  field  strength  can  be  changed  by  moving  a  magnet  toward  or  away  from  the 
coil,  by  moving  the  coil  into  or  out  of  the  magnetic  field,  or  by  rotating  the  coil  relative  to 
the  magnet.  Faraday's  law  is  a  fundamental  relationship  derived  from  Maxwell’s 
equations  and  serves  as  a  succinct  summary  of  the  ways  a  voltage  may  be  generated  by  a 
changing  magnetic  environment.  The  induced  emf  in  a  coil  can  be  expressed  as 

^  A(BA)  AO 

Emf  =  -N  =  -  A -  (20) 

At  At 

where  N  is  the  number  of  turns,  O  or  BA  is  the  magnetic  flux,  B  is  the  external 
magnetic  field  and  A  is  the  area  of  the  coil.  The  induced  emf  involves  the  interaction  of 
charge  with  magnetic  field  [8]. 

b.  Induction  motor  under  V/f  Control 

With  regard  to  rated  speed,  adherence  to  the  V/f  principle  means 
increasing  the  stator  voltage  above  the  rated  value,  which  is  not  permitted.  With 
frequencies  above  rated  speed,  the  applied  voltage  is  maintained  constant  at  the  rated 
level,  which  is  known  as  field  weakening.  Here,  the  intensity  of  the  magnetic  field  in  the 
induction  motor  decreases  with  a  frequency  increase  (unaccompanied  by  a  proportional 
voltage  increase).  The  maximum  permitted  frequency  is  that  which  causes  the  motor  to 
rotate  with  the  maximum  permitted  speed.  At  low  frequencies,  voltage  must  be  boosted 
to  compensate  for  the  effects  of  the  stator  resistance  [7]. 

Using  the  parameters  from  Figure  4,  we  can  express  the  electromagnetic 

torque  as 
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where  is  the  rated  frequency,  co^  is  the  electrical  frequency,  is  the  mutual 

reactance,  s  is  the  slip  (5  =  is  the  rotor  resistance  scaled  by  the  turns 

ratio,  and  is  the  instantaneous  phase  A  current  on  the  stator  side.  The  reactance  X'^  is 

the  sum  of  the  leakage  reactance  on  the  rotor  side  (scaled  by  the  turns  ratio)  X^^  and  the 

mutual  reactance  [2],  Torque  curves  (for  various  values  of  frequency  and  voltage)  and 
applied  voltage  as  a  function  of  frequency  are  displayed  in  Figure  6. 
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Figure  6.  V/f  control  of  induction  motor,  after  [7]. 
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c.  Open  loop  V/f  Control 

We  assume  that  the  three-phase  voltages  at  the  stator  terminals  , 

and  are  given  by 


and 


Van=V,„C0s{2Kf-t), 

(  iTl'^ 

Inf-t-  — 

V  J 


V  -V  cos 

cn  m 


iTzf-t - - 

3 


(22) 

(23) 


(24) 


respectively,  where  the  amplitude  and  the  frequency  /  are  variables.  The  V/f 
control  algorithm  implies  that  there  is  a  relationship  between  the  amplitudes  of  the 
voltage  and  the  frequency.  This  relationship  produces  a  constant  ratio  {K).  We  can  see  a 
top-level  view  of  an  open  loop  V/f  controller  in  Figure  7  [5]. 


Figure  7.  Open  loop  V/f  control,  from  [5]. 
d.  Closed  loop  V/f  Control 

An  improvement  of  open  loop  V/f  control  is  closed-loop  speed  control  by 
slip  regulation  as  shown  in  Figure  8. 
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Figure  8.  Closed-loop  speed  eontrol  utilizing  V/f  and  slip  regulation,  from  [9]. 


The  variables  for  Figure  8  are  defined  in  Table  5. 


Table  5.  Variable  definitions  for  Figure  8. 


Variable 

Variable  Name 

co; 

Rated  speed  command 

ry. 

Speed  feedback  signal 

K 

Slip  speed  command 

Slip  frequency  command 

IM 

Induction  Machine 

Constant 

K, 

Constant 

s 

slip 

K 

Stator  voltage 

A  proportional-integral  (PI)  eontroller  ean  be  employed  to  improve  speed 

aeeuraey  by  regulating  slip  speed  of  the  motor  to  keep  the  motor  speed  at  its  set  value. 

The  major  eomponents  eonsist  of  a  DC  souree,  a  three-phase  inverter  and  an  induetion 

motor  with  load.  The  speed  loop  error  generates  the  slip  eommand  a)*j  through  the  PI 
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controller  and  limiter.  The  slip  is  added  to  the  feedbaek  speed  signal  to  generate  the 
frequeney  command  as  shown.  The  slip  frequency  command  cOg  also  generates  the 
voltage  command  through  a  V/f  function  generator,  which  incorporates  the  low  frequency 
stator  drop  compensation. 

A  step  increase  in  the  slip  frequency  command  produces  a  positive  speed 
error,  and  the  slip  eommand  is  set  to  its  maximum  value.  The  drive  aceelerates  at  the 
permissible  inverter  current,  producing  the  maximum  available  torque  until  the  speed 
error  is  reduced  to  a  small  value.  Finally,  the  drive  settles  to  a  slip  speed  for  whieh  the 
motor  torque  balances  the  load  torque.  A  step  decrease  in  the  slip  frequeney  command 
produces  a  negative  speed  error.  The  slip  speed  is  set  to  the  maximum  negative  value, 
and  the  drive  aecelerates  under  regenerative  braking.  It  eontinues  to  do  so  until  the  speed 
error  is  redueed  to  a  small  value  meanwhile  produeing  the  maximum  available  breaking 
torque.  The  operation  then  shifts  to  motoring,  and  the  drive  settles  at  the  slip  speed  for 
whieh  the  motor  torque  equals  the  load  torque  [9]. 

B,  FIELD-ORIENTED  CONTROL 

In  the  1990s,  the  field  of  eontrolled  electrieal  drives  experienced  rapid  expansion 
due  to  the  advanees  of  semiconductors  in  power  electronics  culminating  in  better 
mieroeleetronic  mieroprocessors  and  digital  signal  proeessors.  These  teehnologieal 
advaneements  enabled  development  of  highly  effective  AC  drive  control  with  lower 
power  dissipation  and  more  aceurate  control  structures.  Drive  eontrols  beeame  more 
accurate  beeause  not  only  were  the  DC  current  and  voltage  controlled  but  the  three-phase 
eurrents  and  voltages  were  managed  by  vector  controls.  The  most  common  form  of 
vector  control  is  FOC.  Three  major  points  establish  this  seheme: 

•  Machine  eurrent  and  voltage  spaee  veetors 

•  Transformation  of  a  three-phase  speed  and  time-dependent  system  into  a 
two  coordinate  time-invariant  system 

•  Effeetive  pulsewidth-modulation  (PWM). 

As  a  result  of  these  factors,  FOC  retains  all  advantages  of  DC  machine  control 
and  frees  itself  from  the  mechanical  commutation  drawbacks.  Furthermore,  FOC  leads  to 
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high  dynamic  performance  in  terms  of  response  times  and  power  conversion  by 
delivering  accurate  steady-state  and  transient  control.  A  basic  FOC  scheme  is  shown  in 
Figure  9  [10]. 


The  variables  for  Figure  9  are  defined  in  Table  6. 


Table  6.  Variable  definitions  for  Figure  9. 


Variable 

Variable  Name 

Reference  rotor  flux 

Estimated  rotor  flux 

r 

Reference  torque 

T 

Estimated  torque 

c 

Reference  d -axis  stator  current 

hd 

Estimated  d -axis  stator  current 

.* 

h. 

Reference  q  -axis  stator  current 

4, 

Estimated  q  -axis  stator  current 

D  -axis  stator  voltage 

Q  -axis  stator  voltage 

4 

Instantaneous  stator  current 

Instantaneous  stator  voltage 

IM 

Induction  Machine 

SVM 

Space  Vector  Modulation 
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In  the  majority  of  classic  AC  drives,  the  generation  of  the  three  sine  waves  is 
based  on  motor  electromechanical  characteristics  and  on  an  equivalent  model  for  the 
motor  in  its  steady-state.  The  control  resembles  that  of  three  separate  phase  system 
controls  rather  than  one  control  of  a  three-phase  system,  presenting  some  major 
drawbacks  including; 

•  The  machine  models  and  characteristics  are  valid  only  in  steady-state. 
This  leads  to  allowing  high  peak  voltage  and  current  transients,  which 
damages  dynamic  performance  and  power  conversion  efficiency. 

•  Considerable  difficulty  in  controlling  the  variables  with  sinusoidal 
references;  PI  controllers  are  incapable  of  performing  sinusoidal 
regulation  without  damaging  the  sinusoidal  reference  and  hysteresis 
controllers  introduce  high  bandwidth  noise  into  the  system. 

•  No  three  phase-system  imbalance  management  and  no  consideration  of  the 
phase  interactions. 

•  Lastly,  the  control  structure  must  be  dedicated  according  to  motor  type 
(synchronous  or  induction)  [10]. 

1.  FOC  Concept 

FOC  consists  of  controlling  the  stator  currents  represented  by  a  vector.  It  is  based 
on  projections  which  transform  a  three-phase  speed  and  time-dependent  system  into  a 
two  coordinate  {d  and  q  coordinates)  time-invariant  system.  These  projections  led  to  a 
structure  similar  to  that  of  DC  machine  control.  FOC  machines  require  two  constants  as 
input  references:  the  torque  component  (aligned  with  the  ^-axis)  and  the  flux  component 
(aligned  with  the  r/-axis).  As  FOC  is  based  on  projections,  the  control  structure  handles 
instantaneous  electrical  quantities  gracefully,  enabling  accurate  control  in  all  operating 
conditions.  Thus,  FOC  solves  classic  scheme  problems  by  reaching  a  constant  reference 
for  both  the  torque  and  flux  components  of  the  stator  current.  In  addition,  maintaining  a 
constant  amplitude  of  the  rotor  flux  results  in  a  linear  relationship  between  torque  and 
torque  component  i^q.  We  can  then  control  the  torque  by  controlling  the  torque 
component  of  the  stator  current  vector,  as  can  be  seen  from 

mccy/R-i^^  (25) 

where  m  is  the  torque  and  is  proportional  to  the  rotor  flux  multiplied  by  the  torque 
component  [11]. 
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2,  Choosing  a  Reference  Frame 

In  general,  there  are  three  possible  seleetions  for  the  flux-linkage  vector:  stator 
flux  linkage  vector,  rotor  flux  linkage  vector,  or  magnetizing  flux  linkage  vector.  In  these 
three  cases  the  instantaneous  electromagnetic  torque  can  be  expressed  as 

K  ~  \¥s  I  Ky  (stator  flux  oriented  control)  (26) 

K  ~  \¥r  I C  (rotor  flux  oriented  control)  (27) 

K  -  r’lm  \¥m  I C  (magnetizing  flux  oriented  control)  (28) 

where  for  linear  magnetic  conditions  and  are  constants  and  and 

\4)rn\  are  the  modulus  of  the  stator,  rotor,  and  magnetizing  flux  linkage  space  vectors, 
respectively.  Moreover,  the  torque-producing  stator  currents  in  the  stator,  rotor,  and 
magnetizing  flux  oriented  reference  frame  are  denoted  by  i|y,  ijy,  and  i^,  respectively. 
The  superscript  denotes  the  reference  frame  used.  These  torque-producing  stator  currents 
take  the  role  of  the  armature  current  in  a  DC  drive. 

Equations  (26)-(28)  can  be  derived  from  a  single  equation,  according  to  which 
electromagnetic  torque  for  an  induction  machine  can  be  expressed  as  the  cross-vector 
product  of  the  stator  flux  linkage  and  current  space  vectors  regardless  of  the  reference 
frame  used.  It  follows  from  Equation  (26)  that  when  the  stator  flux  linkage  modulus  is 
constant  and  the  torque-producing  stator  current  changes  rapidly,  quick  torque  response  is 
obtained.  These  same  considerations  hold  for  Equations  (27)  and  (28).  This  physical 
picture  is  the  basis  for  vector-controlled  drives. 

The  expression  for  the  electromagnetic  torque  contains  the  transformed  stator 
currents,  and  in  a  vector-controlled  drive  the  stator  currents  must  clearly  be  transformed 
into  the  required  reference  frame.  A  common  feature  of  all  vector-controlled  drives  is  that 
the  modulus  and  phase  angle  of  the  AC  excitation  are  controlled.  This  is  the  essence  of 
vector  control.  The  reference  frame  is  aligned  with  the  selected  flux  linkage  space  vector 
so  the  transformation  contains  the  angle  of  the  flux  linkage  space  vector.  This 
implementation  depends  upon  the  flux  angle  and  the  flux  linkage  modulus,  and  a  major 
task  for  implementation  requires  accurate  flux  linkage  estimation  [4]. 
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3,  Defining  the  Space  Vector 

The  three-phase  voltages,  eurrents,  and  fluxes  of  AC  motors  can  be  analyzed  in 
terms  of  complex  space  vectors.  With  regard  to  the  currents,  the  space  vector  are  defined 
as  follows  when  i^,  ib,  and  ic  are  defined  as  the  instantaneous  currents  in  the  stator 
phases.  The  complex  stator  current  ig  is  defined  by 

4  =4+ <^4+ <^"4  (29) 

.2  .4 

where  oc=  and  oc^=  represent  the  spatial  operators.  Shown  in  Figure  10  is  the 
stator  current  complex  space  vector  where  {a,b,c)  are  the  three-phase  system  axes. 


Figure  10.  FOC  of  3-phase  AC  motors,  from  [11]. 


This  current  space  vector  depicts  the  three-phase  sinusoidal  system  that  needs  to  be 
transformed  into  a  two-coordinate  time-invariant  system.  This  transformation  is  a  two- 
step  process: 

•  The  Clarke  transformation,  which  outputs  a  two-coordinate  time-invariant 
system  (a,&,c)=>(oc,  yS)  and 

•  The  Park  transformation,  which  outputs  a  two-coordinate  time-invariant 
system  (oc,  yg)  =>  (J,  g) . 
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a.  Clarke  Transformation 

The  current  space  vector  can  be  reported  in  another  reference  frame  with 
two  orthogonal  axis  called  (oc,  )ff).  Assuming  that  the  axis  a  and  the  axis  oc  are  in  the 
same  direction,  we  have  the  vector  diagram  illustrated  in  Figure  1 1 . 


Figure  1 1 .  Clarke  transformation,  from  [11]. 


The  Clarke  transformation  is  a  projection  that  modifies  the  three-phase 
system  into  the  (oc,  p)  two-dimension  orthogonal  system  and  is  represented  by 

=  4  (30) 

and 


.  _  1  .  2  . 


(31) 


and,  thus,  we  obtain  a  two-coordinate  system  (isoc,  isfi)  that  still  depends  on  time  and 
speed. 
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b.  Park  Transformation 

The  most  important  transformation  in  the  FOC  is  the  Park  transformation. 
This  projeetion  modifies  a  two-phase  orthogonal  system  into  the  {d,q)  rotating 

reference  frame.  If  we  consider  the  (/-axis  aligned  with  the  rotor  flux  for  the  current 
vector,  the  relationship  from  the  two  reference  frames  is  shown  in  Figure  12. 


Knowledge  of  the  rotor  flux  position  6  allows  for  the  flux  and  torque  components 
to  be  determined  by 


hci=h^  cos  6* +  4^  sin^ 


(32) 


and 


sin  6' +  4^  cos  6' 


(33) 


where  and  ig^  are  the  flux  and  torque  components,  respectively.  These  components 

depend  on  the  current  vector  (oc,  )ff)  components  and  on  the  rotor  flux  position.  If  the 
rotor  flux  position  is  known,  then  the  {d,q)  component  is  a  constant.  We  obtain  a  two- 
coordinate  system  (isddsq)  with  the  following  characteristics:  a  two  coordinate  time- 
invariant  system,  and  with  is^  and  igq  direct  control  becomes  possible. 
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c. 


Inverse  Park  Transformation 


The  inverse  Park  transformation  introduces  a  voltage  transformation  that 
modifies  the  voltages  in  the  {d,q)  rotating  reference  frame  in  a  two  phase  orthogonal 
system.  This  transformation  is  represented  by 

^S^ref  =  ^Sdref  COS  0  -  V Sin  0  (34) 


and 


^Sflref  =  ^Sdref  Sin  0  +  COS  0 


(35) 


where  the  outputs  of  this  block  are  the  components  of  the  reference  vector  that  we  call  V^., 
the  voltage  space  vector  to  be  applied  to  the  motor  phases. 


d.  Putting  It  All  Together 

The  basic  FOC  torque  control  scheme  is  summarized  in  Figure  13. 


The  variables  for  Figure  13  are  defined  in  Table  7. 
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Table  7.  Variable  definitions  for  Figure  13. 


ariable 

Variable  Name 

^Sqref 

Flux  reference  stator  current 

^Sdref 

Torque  reference  stator  current 

'sq 

Park  transformed  current 

hd 

Park  transformed  current 

ha 

Clarke  transformed  current 

Clarke  transformed  current 

4 

Phase  A  current 

h 

Phase  B  current 

V 

^  Sqref 

Flux  reference  stator  voltage 

V 

Sdref 

Torque  reference  stator  voltage 

V 

^  Saref 

Inverse  Park  transformed  voltage 

V 

^  Spref 

Inverse  Park  transformed  voltage 

V 

^  DC 

DC  Voltage 

Two  motor  phase  eurrents  are  measured  and  feed  the  Clarke 
transformation  module.  The  outputs  of  this  projeetion  are  designated  isoc  and  isp.  These 
two  components  of  the  current  are  the  inputs  of  the  Park  transformation  that  projects  the 
current  in  the  (d,  q)  rotating  reference  frame.  The  and  components  are  compared 
to  the  references  (the  flux  reference)  and  (the  torque  reference),  respectively. 

This  control  structure  can  be  used  to  control  either  synchronous  or 
induction  machines  by  simply  changing  the  flux  reference  and  obtaining  the  rotor  flux 
position.  In  synchronous  permanent  magnet  motors,  the  rotor  flux  is  fixed  (determined  by 
the  magnets)  and  there  is  no  need  to  create  one.  Hence,  when  controlling  a  permanent 
magnet  synchronous  machine,  should  be  set  to  zero.  However,  an  induction  motor 

needs  a  rotor  flux  to  operate,  and  the  flux  reference  must  not  be  zero.  This  allows  the 
portability  from  asynchronous  to  synchronous  drives  and  solves  a  major  drawback  of 
classic  control  structures.  The  torque  command  is  the  output  of  the  speed  regulator 

when  we  use  a  speed  FOC.  The  outputs  of  the  current  regulators  are  and 
which  are  the  inputs  of  the  inverse  Park  transformation.  The  outputs  of  this  projection  are 
^s«re/  Vspref^  which  arc  the  components  of  the  stator  vector  voltage  in  the  {a,b) 
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stationary  orthogonal  reference  frame.  These  are  the  inputs  of  the  Space  Vector  PWM, 
and  the  outputs  of  this  block  drive  the  inverter.  Remember  both  the  Park  and  inverse  Park 
transformations  need  the  rotor  flux  position  and  obtaining  this  rotor  flux  position  depends 
on  the  AC  machine  type  (synchronous  or  asynchronous  machine).  Rotor  flux  position 
considerations  are  made  in  a  following  subsection  [11]. 

e.  The  Core  of  FOC 

Knowledge  of  the  rotor  flux  position  is  the  core  of  FOC.  Any  error  in  this 
variable  causes  the  rotor  flux  to  be  misaligned  with  the  d  -axis.  In  turn,  and  are 

incorrect  flux  and  torque  components  of  the  stator  current.  The  correct  position  of  the 
rotor  flux,  the  stator  current  and  stator  voltage  space  vector  that  rotates  with  the  {d,q) 
reference  at  synchronous  speed  are  shown  in  Figure  14. 


b 


c 


Figure  14.  Current,  voltage  and  rotor  flux  space  vectors  in  the  d,q  rotating  reference  frame 
and  their  relationship  with  a,b,c  and  a,b  stationary  reference  frames,  from  [11]. 
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f.  Space  Vector  PWM 

The  structure  of  a  typical  three-phase  power  inverter  is  shown  in  Figure 
15,  where  V^,  Vg,  and  V^  are  the  voltages  applied  to  the  wye-connected  motor  windings, 

and  where  is  the  continuous  inverter  input  voltage. 


These  six  switches  must  respect  the  following  conditions; 

•  Three  of  them  must  always  be  on  and  three  must  always  be  off 

•  The  upper  and  lower  switches  of  the  same  leg  are  driven  with  two 
complementary  pulsed  signals  so  that  no  vertical  conduction  is  possible  as 
long  as  there  is  no  overlap  in  the  power  switch  transitions. 

SVPWM  supplies  the  AC  machine  with  the  desired  phase  voltages.  The 

SVPWM  method  of  generating  the  pulsed  signals  fits  the  above  requirements  and 

minimizes  the  harmonic  content.  This  harmonic  content  determines  the  copper  losses  of 

the  machine,  a  major  portion  of  the  machine  losses.  If  these  two  constraints  are  taken  into 

account,  there  are  eight  possible  combinations  for  the  switch  commands.  These  eight 

switch  combinations  determine  eight  phase  voltage  configurations.  These  combinations 

are  depicted  in  Figure  16. 
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Figure  16.  SVPWM,  vectors  and  sectors,  from  [11]. 


These  vectors  divide  the  plane  into  six  sectors.  Depending  on  the  sector 


that  the  voltage  reference  is  in,  we  choose  two  adjacent  vectors.  The  binary 
representation  of  two  adjacent  basic  vectors  differ  in  only  one  bit,  limiting  inverter 
switching  frequency  when  the  switching  pattern  moves  from  one  vector  to  the  next.  The 
two  vectors  are  time  weighted  in  a  sample  period  T  to  produce  the  desired  output 
voltage. 


g.  An  Example  of  SVPWM 


For  example,  if  the  reference  vector  is  in  the  3°  sector  of  Figure  15, 
we  have  the  situation  shown  in  Figure  17. 


V4  (100) 


Figure  17.  Reference  vector  as  a  combination  of  adjacent  vectors,  from  [11]. 
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The  times  and  are  the  times  during  whieh  the  vectors  and  are  applied,  and 
Tq  is  the  time  during  which  the  zero  vectors  are  applied.  When  the  reference  voltage  and 
the  sample  periods  are  known,  the  following  equations 

T  =  T,  +  T,  +  T,  (36) 

K.f=YK+^K  (37) 

make  it  possible  to  determine  the  unknowns  T^,  T^,  and  .  Under  these  constraints  the 

locus  of  the  reference  vector  is  the  inside  of  a  hexagon  whose  vertices  are  formed  by  the 
tips  of  the  eight  vectors.  The  generated  space  vector  PWM  waveforms  are  symmetrical 
with  respect  to  the  middle  of  each  PWM  period.  The  waveforms  in  the  example  presented 
above  are  pictured  in  Figure  18  [11]. 
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Figure  18.  Pattern  of  SVPWM  in  the  3°  sector,  from  [11]. 
Figure  19  illustrates  the  pattern  of  SVPWM  for  each  sector  [11]. 
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Figure  19.  Flexagon  of  SVPWM  pattern,  from  [11]. 


C.  DIRECT  TORQUE  CONTROL 

Direct  torque  control  (DTC)  combines  FOC  theory,  direct  self-control  theory,  and 
advanced  DSP  and  application  specific  integrated  circuit  (ASIC)  technology  to  achieve  a 
sensorless  variable-frequency  drive.  On  the  basis  of  the  errors  between  the  reference  and 
the  estimated  values  of  torque  and  flux,  it  is  possible  to  directly  control  the  inverter  states 
to  reduce  the  torque  and  flux  errors  within  prefixed  band  limits.  Unlike  FOC,  DTC  does 
not  require  any  current  regulator,  coordinate  transformation  or  a  PWM  modulator  (as  a 
consequence  timers  are  not  required).  In  spite  of  its  simplicity,  DTC  allows  for  good 
torque  control  in  both  steady-state  and  transient  operating  conditions.  On  the  other  hand, 
DTC  presents  some  significant  disadvantages  as  follows: 

•  Difficulty  controlling  the  torque  and  flux  at  very  low  speeds. 

•  Produces  high  current  and  torque  ripple. 

•  Possesses  variable  switching  behavior. 

•  Experiences  high  noise  levels  at  low  speed. 

•  Lacks  of  direct  current  control. 
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Research  has  been  carried  out  to  attempt  to  solve  these  issues  faced  by  the  DTC  scheme. 
In  particular,  the  following  solutions  have  been  developed; 

•  Use  of  improved  switching  tables. 

•  Use  of  comparators  with  and  without  hysteresis. 

•  Implementation  of  DTC  schemes  for  constant  switching  frequency 
operation  with  PWM  or  SVM  techniques. 

•  Introduction  of  fuzzy  logic  or  neuro-fuzzy  logic  techniques. 

•  Use  of  sophisticated  flux  estimators  to  improve  low  speed  behavior. 

Although  all  of  these  approaches  improve  the  performance  of  DTC,  they  also  lead 

to  more  complex  schemes  detracting  from  DTC  simplicity.  The  DTC  technique  is 
intrinsically  sensorless,  so  only  the  control  schemes  which  meet  this  requirement  (i.e.,  no 
current  or  PI  regulators,  no  coordinate  transformations,  and  no  PWM  signals  generator) 
are  considered  true  DTC  schemes  [10]. 

1.  DTC  of  a  VSI-Fed  Induction  Motor 

It  is  possible  to  directly  control  the  stator  flux  linkage  (or  rotor  or  magnetizing 
flux  linkages)  and  the  electromagnetic  torque  by  the  selection  of  optimum  inverter 
switching  modes.  The  selection  is  made  to  restrict  the  flux  and  torque  errors  within 
respective  flux  and  torque  hysteresis  bands,  to  obtain  fast  torque  response,  low  inverter 
switching  frequency,  and  low  harmonic  losses.  DTC  allows  very  fast  torque  responses 
and  flexible  control  of  an  induction  machine  [4]. 

2.  General,  Mathematical,  and  Physical  Fundamentals  of  Producing 
Fast  Torque  Response 

a.  Stator  Flux  Linkage,  Rotor  Flux  Linkage,  and  Stator  Current 
Space  Vectors 

In  a  balanced  three-phase  induction  machine,  torque  is  proportional  to  the 
cross-product  of  stator  flux  linkage  and  stator  current  space  vectors,  as  given  by 

(38) 
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where  is  the  stator  flux  linkage  space  vector  and  4  is  the  stator  current  space  vector. 
In  Equation  (38),  both  space  vectors  are  expressed  in  the  stationary  reference  frames.  By 
considering  that  |  exp(j  ),  where  is  the  angle  of  the  stator  flux  linkage  space 

vector  with  respect  to  the  direct  axis  of  the  stator  reference  frame  and  4  “  |4 1  )■>  we 

can  put  Equation  (38)  into  the  form 

3  —  3  — 

4  =  2^I^J|4|sin(«,-A)  =  -f’|^J|4|sin«  (39) 

where  a  =  -  pps  the  angle  between  the  stator  flux  linkage  and  stator  current  space 

vector.  The  relationship  between  the  stator  flux  linkage  and  stator  current  space  vectors  is 
shown  in  Eigure  20. 
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Eigure  20.  Stator  flux  linkage  and  stator  current  space  vectors,  from  [4]. 

If  the  stator-flux  linkage  space  vector  is  kept  constant  and  p^  is  changed 
quickly,  the  electromechanical  torque  can  also  be  rapidly  changed  {Atjdt  «  Apjdt). 
By  forcing  the  largest  d p^  jdt ,  the  fastest  torque  response  time  is  obtained. 

If  stator  voltages  are  imposed  on  the  motor  which  keep  the  stator  flux 


constant  but  which  quickly  rotate  the  stator  flux  linkage  space  vector  into  the  required 
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position  demanded  by  the  torque,  then  fast  torque  eontrol  is  aehieved.  The  parameter 
ean  be  quiekly  ehanged  by  eontrolling  the  stator  flux  linkage  spaee  veetor  using  the 
appropriate  stator  voltages  (generated  by  the  VSI-fed  induction  machine).  By  using  the 
appropriate  stator  voltages,  we  achieve  direct  stator  flux  and  electromagnetic  torque 
control. 


An  extremely  clear  picture  of  this  process  is  painted  if  we  consider  another 
form  of  Equation  (39),  leading  to  the  same  results  mentioned  above.  By  considering  that 
(p^  =  +LJ^  and  (p\.  =  +  where  the  primed  rotor  quantities  are  expressed  in 

the  stationary  reference  frame,  we  see  that  z)  =  Equation  (39) 

takes  on  the  form 


3  r  ,  3  E 

t  =—P  ,  m  xm  =—P—p^ 
"  2  L.L  "  *  2  L.L 


K||^.|sin(p^-pJ  =  |E^ 


(p\  (p^  siny  (40) 


where  y  is  the  angle  between  the  stator  and  rotor  flux  linkage  space  vectors  and  is  the 

angle  of  the  rotor  flux  linkage  space  vector  with  respect  to  the  real  axis  of  the  stationary 
reference  frame  as  shown  in  Eigure  21. 


Eigure  21.  Stator  flux  linkage,  rotor  flux  linkage,  and  stator  current  space  vectors,  from  [4]. 
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b.  Achieving  Rapid  Torque  Control 

The  rotor  time  eonstant  is  large;  therefore,  the  rotor  flux  linkage  changes 
slowly  compared  to  the  stator.  It  can  be  assumed  to  remain  constant.  This  follows  from 
the  rotor  voltage  equation  if  the  stator  flux  linkage  is  assumed  to  be  constant.  If  both  are 
assumed  to  be  constant,  the  torque  can  be  rapidly  changed  by  changing  y  in  the  required 
direction  (determined  by  the  torque  command).  This  is  the  essence  of  DTC.  By  choosing 
the  appropriate  stator  voltage  space  vector,  the  angle  y  can  be  managed.  Even  if  the 

stator  flux  linkage  is  not  held  constant,  both  y  and  (p^  can  be  controlled.  In  FOC  the 

stator  currents  are  the  control  quantities,  while  in  DTC  the  stator  flux  linkages  are 
controlled.  Equation  (41)  relates  to  that  of  a  synchronous  machine  where  the  torque  is 
controlled  by  the  load  angle  between  the  stator  and  rotor  flux  linkages. 

In  summary,  by  rapidly  changing  the  position  or  speed  of  the  stator  flux 
linkage  space  vector,  rapid  torque  control  can  be  achieved.  Both  the  modulus  and  the 
angle  of  the  stator  flux  linkage  space  vector  can  be  changed  by  the  stator  voltages. 
Correct  step-by-step  selection  of  the  stator  voltage  vector  permits  the  stator  flux  to  be 
changed  in  the  required  way.  Decoupled  control  of  the  torque  and  stator  flux  is  achieved 
by  acting  on  the  radial  and  tangential  components  of  the  stator  flux  linkage  space  vector. 
These  two  components  are  directly  proportional  to  the  components  of  the  stator  voltage 
space  vector  in  the  same  directions;  therefore,  they  can  be  controlled  by  the  appropriate 
inverter  switchings.  The  angle  y ,  or  the  relative  position  of  the  stator  and  rotor  flux 
linkage  space  vectors,  determines  torque.  By  assuming  a  slow  motion  of  the  rotor,  if  a 
stator  voltage  is  applied  which  causes  a  quick  movement  of  the  stator  flux  linkage  away 
from  that  of  the  rotor  flux  linkage,  the  torque  increases  because  y  increases.  On  the  other 
hand,  if  a  zero-voltage  space  vector  is  applied  which  nearly  stops  the  rotation  of  the  stator 
flux  linkage,  then  torque  decreases  because  the  rotor  flux  linkage  is  still  rotating  and  y 
decreases.  If  a  zero-voltage  vector  is  applied  long  enough  (where  stator  flux  linkage 
almost  does  not  move),  the  rotor  flux  linkage  overtakes  that  of  the  stator  and  the  sign  of 
y  changes,  and  torque  changes  direction.  For  these  reasons,  the  duration  that  the  zero- 
voltage  space  vectors  are  applied  have  a  direct  effect  on  torque  oscillations  [4]. 
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c.  PWM  VSI  Inverter  Switching  States  and  Space  Vectors 
A  six-pulse  PWM  VSI-fed  inverter  is  shown  in  Figure  22. 


This  inverter  produces  six  non-zero  active  voltage-switching  space  vectors 
(Wj,  iTj,  and  two  zero  space  vectors  (w^,  u^).  These  eight  switching  states  are 

shown  in  Figure  23,  and  their  corresponding  switching  vectors  are  shown  in  Figure  24. 
The  six  active  inverter-switching  vectors  can  be  expressed  as 

(42) 

where  is  the  DC  link  voltage  and  A:  =1,2, ...6.  For  k=l  and  8,  %  =0holds  for  the 
two  zero-switching  states  where  the  stator  windings  are  short  circuited,  and  =  0  . 
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u,  (100)  U2(110)  iiaCOlO)  m4(011) 


a  b  c  a  b  c  a  b  c  a  b  c 


Us  (001)  «6(101)  «7(111)  58(000) 


a  b  c  a  b  c  a  b  c  a  b  c 


Figure  23.  Eight  switching  states  of  the  PWM  VSI  inverter,  from  [4], 


Figure  24.  Eight  corresponding  switching  space  vectors,  from  [4]. 


d.  Controlling  Torque  and  Flux  through  Hysteresis 

During  every  sampling  period,  we  select  switching  vectors  to  keep  the 
stator  flux  linkage  errors  within  the  required  tolerance  band  and  the  torque  error  within 
the  hysteresis  band.  The  widths  of  these  bands  are  2Atg  and  2A^^ ,  respectively.  The 
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upper  limit  is  one  A  above  the  referenee,  and  the  lower  limit  is  one  A  below.  In  terms  of 
low  order  harmonics,  the  flux  hysteresis  band  mainly  affects  the  stator  current  distortion. 
The  torque  hysteresis  band  affects  the  switching  frequency  and,  thus,  the  switching 
losses.  This  scheme  requires  the  use  of  flux  linkage  and  torque  estimators.  Stator  flux 
linkage  components  can  be  obtained  by  integrating  appropriately  monitored  terminal 
voltages.  On  the  other  hand,  it  is  unnecessary  to  monitor  stator  voltages  because  they  can 
be  reconstructed  by  using  the  inverter  switching  modes  and  monitored  DC  link  voltage. 

The  goal  is  to  keep  the  modulus  of  the  stator  flux  linkage  space  vector  | 
within  the  hysteresis  band,  denoted  by  the  two  circles  whose  width  is  2A^^  as  shown  in 
Figure  25. 
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Figure  25.  Control  of  the  stator  flux  linkage  space  vector:  stator  flux  linkage  space  vector 

locus  and  inverter  switching  errors,  from  [4]. 
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Due  to  the  six-step  inverter,  the  loeus  of  the  flux  linkage  spaee  veetor  is  divided  into  the 
six  sectors  shown  in  Figure  24.  If  the  stator  flux  linkage  space  vector  is  initially  at  point 
(sector  one)  and  rotating  counterclockwise,  then  the  stator  flux  linkage  is  at  the  upper 

limit  and  must  be  reduced.  As  shown  in  Figure  25,  applying  achieves  this  reduction. 
The  stator  flux  linkage  space  vector  moves  from  to  (sector  2)  and,  once  again, 
reaches  the  upper  limit.  Once  again  (when  rotating  counterclockwise),  applying  the 
switching  vector  reduces  the  stator  flux  linkage  appropriately.  Now,  moves  from 

Py  io  P2  (still  in  sector  2).  If  at  point  Py ,  the  rotation  of  the  stator  flux  linkage  space 

vector  must  be  stopped,  then  a  zero-switching  vector,  either  Uj  (111)  or  (000),  has  to 

be  applied.  Since  the  last  switching  vector  applied  was  (010),  then  is  used  because 

only  the  switching  of  the  second  switch  from  1  to  0  is  required.  The  switching  vector  that 
requires  the  minimum  amount  of  switching  should  always  be  chosen  [4]. 

As  discussed  above,  stopping  the  rotation  of  the  stator  flux  linkage  space 
vector  corresponds  to  the  case  when  the  electromagnetic  torque  does  not  have  to  be 
changed.  When  it  has  to  be  changed  in  either  the  clockwise  or  counterclockwise 
direction,  then  the  stator  flux  linkage  space  vector  has  to  be  rotated  in  the  appropriate 
direction.  When  the  stator  flux  rotates  counterclockwise  and  if  an  increase  in  torque  is 
required  (i.e.,  the  stator  flux  linkage  space  vector  is  at  point  Py  where  the  flux  linkage  is 
at  its  upper  limit),  then  the  increase  in  torque  is  achieved  by  applying  switching  vector 
if, .  If  a  decrease  in  torque  is  required  but  the  flux  linkage  is  at  the  lower  limit  and  must 

be  increased,  then  this  is  achieved  by  applying  switching  vector  iq  because  this  moves 
the  stator  flux  linkage  space  vector  in  the  clockwise  direction  (or  in  the  direction  of 
negative  torque)  and  also  increases  the  flux  linkage. 

e.  Choosing  the  Optimum  Voltage  Switching  Vector 

The  positions  of  the  various  stator  flux  linkage  vector,  if  the  stator  flux 
linkage  space  vector  is  in  one  of  the  six  sectors  are  shown  in  Figure  26. 
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V|4  in  Sector  6 

Figure  26.  Position  of  various  stator  flux  linkage  spaee  veetors  and  seleetion  of  the  optimum 
switehing  veetors.  FI;  flux  inerease;  FD;  flux  deerease;  TI;  torque  inerease;  TD; 

torque  deerease,  from  [4]. 


In  general,  if  an  inerease  of  the  torque  is  required,  then  the  torque  is 

eontrolled  by  applying  voltage  veetors  that  advanee  the  flux  linkage  spaee  veetor  in  the 

direetion  of  rotation.  If  a  deerease  is  required,  voltage  veetors  are  applied  whieh  oppose 
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the  direction  of  torque.  If  zero  torque  is  required,  the  zero-switching  vector  is  applied  that 
minimizes  inverter  switching.  The  angle  of  the  stator  voltage  space  vector  is  indirectly 
controlled  through  the  flux  vector  modulus  and  torque,  and  increasing  torque  results  in  an 
increased  angle.  The  torque  demand  is  reduced  to  a  choice  of  increase,  decrease,  or  zero 
while  the  stator  flux  linkage  vector  modulus  can  either  be  increased  or  decreased.  These 
results  are  tabulated  in  the  optimum  switching  vector  selection  table  as  shown  in  Table  8 

[4]. 


Table  8.  Optimum  voltage  switching  vector  look-up  table,  from  [4]. 


cl(/r 

di. 

a(l) 

sector  1 

a(2) 
sector  2 

a(3) 

sector  3 

a(4) 

sector  4 

ot(5) 
sector  5 

a(6) 

sector  6 

1 

“2 

«3 

“5 

“6 

Ml 

1 

0 

-1 

«8 

“7 

“8 

M7 

Ms 

«6 

“2 

“4 

M$ 

1 

“3 

“5 

“6 

Ml 

M2 

0 

0 

“8 

“7 

Ws 

Mfl 

M7 

«5 

«6 

Ml 

“2 

M3 

M4 

Active  switching  vectors:  u,(I00);  Uj(l  10);  UjlOlO);  U4(0l  1);  iJj(OOl); 
Zero  switching  vectors:  U7(l  1 1);  UglOOO). 


The  optimum  selection  of  the  switching  vectors  for  all  the  possible  stator 
flux  linkage  space  vector  positions  and  the  desired  control  inputs  (which  are  the  reference 
values  of  the  stator  flux  linkage  modulus  and  the  electromagnetic  torque,  respectively) 
are  given  in  Table  8.  The  optimum  switching  look-up  table  requires  knowledge  of  the 
position  of  the  stator  flux-linkage  space  vector  because  it  must  be  known  in  which  sector 
it  resides.  For  this  purpose,  the  angles  a(V) ,  «(2),...,«(6)  shown  in  Figure  24  must  be 
known.  We  can  eliminate  the  need  for  trigonometric  functions  because  it  is  not  the 
accurate  position  of  the  stator  flux  linkage  space  vector  that  has  to  be  known,  but  only  the 
sector  in  which  it  resides.  This  information  can  be  obtained  by  considering  the  signs  of 
the  various  stator  flux  linkage  components,  allowing  a  simple  implementation  requiring 
only  the  use  of  comparators.  The  results  of  these  signs  are  summarized  in  Table  9  [4]. 
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Table  9.  Selection  of  the  stator  flux  linkage  space  vector  sector,  from  [4], 


Sectors 

Signs  of  flux  linkages  ^ 

«(I)  a(2)  2(3) 

sector  1  sector  2  sector  3 

«(4) 
sector  4 

*(5) 
sector  5 

2(6) 

sector  6 

Sign  of  1^,0 

+  +  - 

- 

- 

+ 

Sign  of  4i,q 

Sign  of 

{nu;-+)  +  + 

(nu;  +  -) 

" 

=  sign  of 

-  +  + 

+ 

— 

— 

nu  -  not  useful 


The  relationship  of  the  space  vector  (p^  to  the  various  stator  flux  linkage  components  is 
shown  in  Figure  27. 
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(sA) 


Figure  27.  Relationship  of  the  space  vector  (p^  to  the  stator  flux  linkage  components  (p^,^  , 

(p^Q  ,  and  (p^B  ,  from  [4]. 


It  should  be  noted  that  at  very  low  speeds,  such  as  when  the  machine  is 
started,  flux  control  can  be  lost  during  application  of  the  switching  vectors  shown  in 
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Table  5.  Instead  of  having  a  circular  path  corresponding  to  constant  reference  flux,  the 
space  vector  locus  of  the  stator  flux  linkage  space  vector  is  a  six-sided  symmetrical  locus, 
where  during  1/6  of  the  cycle  the  modulus  of  the  stator  flux  linkage  space  vector  changes. 
These  problems  are  related  to  the  inappropriate  use  of  switching  voltage  vectors  in  the 
low  speed  region. 

f.  A  Basic  DTC  Scheme 

As  can  be  seen  in  Figure  28,  the  error  between  the  estimated  torque  T  and 

* 

the  reference  torque  T  is  the  input  of  a  three-level  hysteresis  comparator  which  is  shown 
in  Figure  29.  The  error  between  the  estimated  stator  flux  magnitude  (Pg  and  the  reference 
stator  flux  magnitude  (p^,  is  the  input  of  a  two-level  hysteresis  comparator,  which  is 
shown  in  Figure  30  [10]. 


The  variables  in  Figure  28  are  defined  in  Table  10. 
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Table  10.  Variable  definitions  for  Figure  28. 


Variable 

Variable  Name 

>l< 

(Ps 

Refemce  stator  flux  magnitude 

(Ps 

Estimated  stator  flux  magnitude 

T* 

Reference  torque 

T 

Estimated  torque 

Output  of  flux  hysteresis  comparator 

Ct 

Output  of  torque  hysteresis  comparator 

IM 

Induction  Machine 

Figure  29.  Torque  hysteresis  comparator,  from  [10]. 
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Figure  30.  Flux  hysteresis  comparator,  from  [10]. 


The  selection  of  the  appropriate  voltage  vector  is  based  on  the  switching 
table  given  in  Table  11.  Knowledge  of  the  optimum  switching  look-up  table  requires 
knowledge  of  the  stator  flux  linkage  space  vector  since  it  must  be  known  in  which  sector 
the  vector  is  [4]. 
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Table  11.  Basic  switching  table  for  DTC,  from  [10]. 


Sector 

1 

2 

3 

■■ 

5 

6 

- 1 

S  =  -l 

\  Ct  - 

^5 

h 

o 

11 

K 

V7 

Cj  =  +1 

V, 

V, 

V, 

Vs 

■e 

II 

+ 

Cj  — 

V, 

Ve 

K 

h 

o 

II 

K 

^7 

K 

Vi 

m 

V, 

I  Cj-  —  +1 

b 

lAj 

m 

K 

The  input  quantities  are  the  stator  flux  sector  and  the  outputs  of  the  two 
hysteresis  comparators.  If  the  stator  flux  vector  is  lying  in  sector  one  of  the  {d,q)  plane, 
then  the  voltage  vectors  used  by  DTC  technique  are  shown  in  Figure  31. 


d-oxis 


Figure  3 1 .  Voltage  vectors  utilized  in  basic  DTC  scheme  when  stator  flux  in  sector  1 , 

from  [10]. 
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This  simple  approach  allows  a  quick  torque  response  to  be  achieved; 
however,  the  steady-state  performance  is  characterized  by  undesired  ripple  in  current, 
flux  and  torque.  This  is  largely  due  to  the  absence  of  information  about  torque  and  rotor 
speed  values  in  the  voltage  vector  selection  algorithm  [10]. 

g.  Limitations  of  DTC 

(1)  Torque  and  Flux  Ripple.  In  DTC  drives  there  are  torque 
and  flux  ripples  because  none  of  the  inverter  switching  vectors  are  able  to  generate  the 
exact  stator  voltage  required  to  produce  the  desired  changes  in  torque  and  stator  flux 
linkage.  Various  schemes  can  be  used  to  reduce  these  ripples.  Although  increased 
switching  frequency  is  desirable  (because  it  reduces  harmonic  content  of  stator  currents 
and  leads  to  reduced  torque  harmonics),  it  significantly  increases  switching  losses, 
leading  to  reduced  efficiency  and  added  stress  on  the  semiconductor  devices  of  the 
inverter.  In  the  case  of  high  switching  frequency,  a  faster  processor  is  required,  increasing 
cost.  In  addition,  changed  inverter  topology  can  be  used  to  limit  ripples. 

(2)  Problems  at  Low  Speed.  When  a  DTC  drive  is  operated  in 
the  zero-speed  region,  problems  occur.  During  magnetization,  the  stator  flux  comparator 
selects  only  non-zero  vectors.  At  that  time,  the  output  of  the  torque  comparator  takes  one 
state,  and  the  inverter  cannot  apply  zero  vectors  to  the  motor.  The  introduction  of  an 
additional  carrier  signal  to  the  input  of  the  torque  comparator  is  one  solution  to  this 
problem.  This  carrier  signal  is  only  applied  in  the  zero-speed  region  and  forces  zero 
vectors  improving  stator  flux  linkage  and  stator  current  waveforms.  This  ensures  robust 
start  and  operation  in  the  zero-speed  region. 

(3)  Duty  Ratio  Control.  In  DTC  motor  drives,  a  voltage  vector 
is  applied  for  the  entire  switching  period,  which  causes  stator  current  and  torque  to 
increase  for  that  entire  period.  For  small  errors,  torque  exceeds  its  reference  value  early 
and  continues  to  increase,  causing  a  high  torque  ripple.  Then,  zero-vectors  are  applied  to 
reduce  the  torque  back  to  its  reference  value.  A  solution  can  be  obtained  to  reduce  the 
ripples  by  applying  the  inverter  switching  vector  for  only  a  part  of  the  switching  period 
and  a  zero-vector  for  the  remainder  of  the  period,  thus,  a  non-zero  vector  is  used  to 

increase  torque  to  its  reference  value  and,  when  it  reaches  that  value,  a  zero-vector  is 
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applied.  During  the  applieation  of  a  zero-veetor,  there  is  no  voltage  foreed  on  the 
maehine  and  torque  is  nearly  eonstant.  The  average  input  voltage  during  applieation  is 
5U^  (5  is  the  duty  ratio  from  0  to  1).  So  any  voltage  between  0  and  can  be  applied. 

This  inereases  the  number  of  ehoiees  for  the  voltage  veetor,  whieh  was  previously  limited 
to  the  six  switehing  veetors  [4]. 
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IV.  MEASUREMENTS  AND  SIMULATIONS 


A.  APPROACH 

It  is  possible  to  model  virtually  every  aspeet  of  an  induction  machine,  but  this 
model  serves  only  as  a  means  to  reduce  time  spent  in  the  laboratory,  cost  and  risk  to 
personnel.  A  useful  model  must  be  able  to  predict  laboratory  results  and  circuit  dynamics. 
Only  through  experimental  validation  is  it  apparent  that  a  simulated  model  predicts  the 
behavior  of  an  implemented  hardware  prototype  accurately.  Upon  verifying  the  V/f 
control  simulation  functioned  properly,  we  used  this  model  was  used  to  simulate  FOC 
and  DTC.  Ultimately,  three  types  of  torque  control  for  a  VSI-fed  induction  machine  were 
compared. 

The  experimental  data  presented  here  was  captured  through  the  use  of  Chipscope, 
an  embedded,  software  based  logic  analyzer.  This  program  allowed  monitoring  of  all  of 
the  signals  in  the  design  and  provided  a  convenient  software  based  interface  for 
controlling  the  integrated  logic  analyzer,  including  setting  triggering  options  and  viewing 
the  waveforms.  The  program  also  provided  a  fast,  easy,  and  interactive  setup  and  simple 
debugging  of  input/output  channels  in  high-speed  FPGA  designs.  Up  to  4000  data  points 
were  captured  from  four  separate  channels.  The  data  collected  included  the  flux  (CHI), 
Phase  B  current  in  amps  (CH2),  commanded  voltage  in  volts  (CHS),  and  the  speed  in  rpm 
(CH4).  These  signals  were  captured  at  the  speed  of  operation  and  brought  out  through  the 
programming  interface.  Captured  signals  were  displayed  and  analyzed  using  the 
Chipscope  Pro  Analyzer  tool. 

B,  V/F  CONTROL 

To  show  that  the  simulation  accurately  predicted  the  behavior  of  the  machine, 
simulated  results  were  compared  with  measured  data.  It  can  be  seen  in  Figure  32  that  the 
speed  response  of  the  measured  and  simulated  results  match  closely. 
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Figure  32.  Measured  and  simulated  speed  response  of  an  induction  machine  with  V/f 

control. 

Although  only  a  certain  number  of  samples  were  gathered  (due  to  limitations  on 
available  lab  equipment),  the  data  that  was  collected  sufficiently  demonstrates  that  the 
simulation  accurately  predicted  the  behavior  of  a  hardware  prototype. 

The  primary  objective  was  to  illustrate  the  operation  of  both  the  physical  machine 
and  the  simulation  in  transient  and  steady-state  operating  conditions.  This  was  achieved 
by  creating  three  distinct  step  events  and  recording  the  results.  The  machine  was 
originally  set  to  80  percent  of  its  rated  speed  (approximately  1440  rpm).  The  first  step 
event  occurs  at  0.5  seconds,  where  the  speed  is  increased  from  80  percent  of  the  rated 
value  to  the  rated  value  itself  (1800  rpm).  Although  the  maximum  rated  value  was  never 
reached  (due  to  slip  and  friction),  neither  the  model  nor  the  machine  itself  reached  this 
speed.  The  second  step  event  occurs  at  two  seconds,  where  the  speed  was  decreased  from 
the  rated  value  back  to  the  original  80  percent  level.  Finally,  the  third  step  event  occurs  at 
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three  seconds,  and  the  speed  was  decreased  from  80  percent  to  60  percent  (approximately 
1080  rpm)  of  the  rated  value.  Throughout  these  step  events,  the  measured  and  simulated 
data  align  and  demonstrate  similar  response  times. 

The  Phase  B  current  during  these  same  step  events  was  also  monitored.  This 
current  was  filtered  at  300  Hz,  and  the  step  response  of  the  step  event  occurring  at  0.5 
seconds  (from  80  percent  to  100  percent  of  the  rated  speed)  is  shown  in  Figure  33. 
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Figure  33.  Phase  B  current  filtered  at  300  Hz  of  an  induction  machine  under  V/f 

control. 


At  steady-state,  the  measured  and  simulated  currents  are  in  phase.  As  this  model 
was  developed  using  steady-state  operating  conditions,  the  two  did  not  obtain  the  same 
transient  response.  This  observation  was  made  for  the  other  two  step  events  as  well.  With 
the  measured  data,  we  see  a  current  spike  of  over  6.0  A,  whereas  the  maximum  simulated 
current  spike  is  only  about  2.3  A.  In  addition,  during  the  transient  response,  the  currents 
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fell  out  of  phase,  and  the  simulated  eurrent  led  the  measured  eurrent.  This  error  was 
predieted  to  oeeur  beeause  the  steady-state  equivalent  eireuit  model  does  not  aeeurately 
model  the  dynamie  response;  however,  the  error  is  small.  Upon  settling  back  to  steady- 
state,  the  behavior  of  the  machine  fell  right  in  line  with  what  the  simulation  predicted. 
Based  on  the  speed  and  current  response  to  the  three  created  step  events,  the  model 
accurately  predicted  the  behavior  of  an  induction  machine  utilizing  V/f  control.  This 
verification  allowed  the  model  to  be  used  for  simulating  FOC  and  DTC. 

C.  V/F  CONTROL  VERSUS  FOC  AND  DTC 

Time  and  money  can  be  saved  through  the  use  of  simulation  because  it  is 
unnecessary  to  physically  implement  all  of  the  control  schemes  covered  in  this  thesis. 
The  developed  Simulink  model  has  proven  to  successfully  demonstrate  the  behavior  of  an 
induction  machine  utilizing  V/f  control  and  can  be  built  upon  to  simulate  FOC  and  DTC 
as  well. 


1.  V/f  Control  Versus  FOC 

The  simulation  was  modified  to  incorporate  FOC  capability.  After  these 
modifications  were  made,  this  control  scheme  was  compared  to  V/f  control  by 
introducing  the  same  step  events  used  earlier.  This  comparison  is  shown  in  Figure  34. 
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Figure  34.  Comparison  of  simulated  V/f  eontrol  and  FOC. 

FOC  possesses  a  mueh  better  dynamie  response  time  (it  reaehes  its  maximum 
speed  mueh  quieker  than  V/f  eontrol).  Moreover,  the  full  rated  speed  of  the  induetion 
maehine  is  aeeomplished  (1800  rpm).  These  two  facts  prove  that  FOC  possesses  distinct 
advantages  over  V/f  control. 

2.  V/f  Versus  DTC 

Next,  DTC  capability  was  developed  and  incorporated  into  the  Simulink  model. 
In  the  same  way  that  FOC  was  compared  to  V/f  control,  so  was  DTC.  This  comparison  is 
shown  in  Figure  35. 
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Figure  35.  Comparison  of  simulated  V/f  control  and  DTC. 


DTC  shows  a  distinct  advantage  over  V/f  control  for  increasing  speeds  but  is 
slightly  sluggish  for  speed  decreases.  Much  like  FOC,  DTC  allows  the  machine  to  reach 
the  full  rated  speed.  Both  the  FOC  and  DTC  simulations  illustrate  clear  benefits  over  V/f 
control. 

D.  FOC  Versus  DTC 

FOC  induction  machines  were  introduced  nearly  40  years  ago  and  today  are  an 
industrial  reality  and  widely  available.  DTC  was  introduced  shortly  thereafter  and  was 
characterized  by  simplicity,  good  performance  and  robustness.  The  use  of  DTC  made 
dynamic  control  of  torque  possible  without  any  mechanical  transducers  on  machine 
shafts.  Often  DTC  is  referred  to  as  a  sensorless  type  control  technique  [1].  It  was  the  aim 
of  this  thesis  to  fairly  compare  these  two  techniques  in  all  modes  of  operation.  First,  DTC 
and  FOC  speed  responses  are  compared  and  are  displayed  in  Figure  36. 


60 


Time  (seconds) 

Figure  36.  Comparison  of  simulated  speed  response  of  DTC  and  FOC. 

Both  eontrol  schemes  reach  the  rated  speed  and  do  so  relatively  quickly. 
However,  FOC  performs  considerably  better  than  DTC  upon  both  decreases  from  100 
percent  to  80  percent  and  from  80  percent  to  60  percent.  It  is  also  important  to  examine 
the  Phase  B  currents  of  both  control  strategies  (still  filtered  at  300  Hz)  for  the  speed 
profiles  shown  in  Figure  36.  They  are  displayed  in  Figure  37. 
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Figure  37.  Phase  B  current  for  DTC  and  FOC. 

DTC’s  peak-to-peak  current  during  steady-state  operation  is  much  lower  than  the 
peak-to-peak  current  of  FOC.  This  lower  current  leads  to  lower  electrical  losses  of  the 
machine.  Although  FOC  constantly  controls  the  current  through  the  use  of  current 
sensors,  DTC  does  not  have  this  same  control  and  experiences  high  current  peaks  during 
transient  operation.  As  discussed  earlier,  DTC  experiences  high  current  ripples  because 
the  optimum  switching  vectors  cannot  generate  the  exact  stator  voltage  required  to 
produce  the  desired  changes  in  torque  and  stator  flux  linkage.  Based  on  these  results, 
each  control  scheme  possesses  advantages,  and  the  application  of  a  machine  determines 
which  one  to  employ.  If  the  machine  is  going  to  be  used  mostly  in  steady-state  operation, 
such  as  a  water  pump,  DTC  is  the  better  choice.  If  the  machine  changes  speeds 
frequently,  such  as  an  electric  drive  for  an  electric  car,  then  the  preferred  choice  is  FOC. 
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Next,  the  reference  torque  of  both  DTC  and  FOC  are  examined  and  shown  in 
Figure  38. 


Time  (seconds) 


Figure  38.  Reference  torque  for  DTC  and  FOC. 


When  considering  the  reference  torque,  it  should  be  noted  that  the  absolute  value 
of  the  reference  torque  was  set  to  two  Newton  meters  as  determined  by  a  saturation  block 
in  the  controller  of  this  model.  It  is  also  important  to  understand  that  the  signed  values  of 
torque  imply  the  demanded  direction  of  torque  where  positive  values  imply  rotation  in 
the  counterclockwise  direction,  while  negative  values  imply  rotation  in  the  clockwise 
direction.  During  the  first  step  event  (at  0.5  seconds),  we  see  that  both  control  schemes 
demand  the  maximum  available  torque  in  the  counterclockwise  direction  in  an  attempt  to 
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speed  the  machine  from  80  percent  of  the  rated  value  to  the  rated  value  itself.  Although 
both  schemes  call  for  torque  in  similar  ways,  FOC  is  noticeably  more  responsive  upon 
speed  decreases. 

While  the  reference  torque  is  the  torque  the  machine  calls  for,  the  electrical  torque 
is  what  the  machine  actually  receives.  The  electrical  torque  is  displayed  in  Figure  39. 


Figure  39.  Electrical  torque  for  DTC  and  FOC. 

The  first  observance  from  the  electrical  torque  is  the  different  shapes  of  both. 
FOC  has  a  pulsed  response,  while  DTC  has  a  stepped  response.  FOC  has  a  slightly 
slower  response  time  than  DTC,  and  the  torque  for  DTC  is  controlled  through  a 
hysteresis  comparator  and,  therefore,  is  bouncing  between  the  upper  and  lower  limits  of 
this  comparator. 

Finally,  we  need  to  consider  how  the  flux  behaves  under  these  two  control 
schemes.  The  stator  flux  of  DTC  and  the  stator  flux  of  FOC  can  be  seen  in  Figure  40. 
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Figure  40.  Simulated  flux  for  DTC  and  FOC. 

The  results  in  Figure  40  were  produeed  by  eapturing  a  few  eyeles  of  the  eleetrieal 
flux  during  steady-state  operation  at  the  rated  speed.  These  eyeles  were  eaptured  from 
1.45  to  1.55  s.  FOC  maintains  eonstant  flux  throughout  the  duration  of  this  eyele 
(produeing  a  cirele),  while  DTC  does  not.  As  discussed  earlier,  DTC  attempts  to  keep  the 
flux  within  a  preset  hysteresis  band  with  upper  and  lower  limits,  and  the  uneven  response 
is  a  result  of  bouncing  back  and  forth  between  these  two  limits.  With  FOC  the  torque  and 
flux  producing  current  components  are  decoupled,  and  the  transient  response 
characteristics  emulate  those  of  a  separately  excited  DC  machine.  Observing  the  dynamic 
response  of  flux  is  also  beneficial.  This  is  shown  in  Figure  41. 
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Figure  41 .  Electrical  flux  for  DTC  and  FOC. 

Again,  DTC  maintains  the  flux  within  the  hysteresis  band.  DTC’s  electrical 
torque  drops  to  zero  and  back  again  because  the  direction  of  the  flux  linkage  space  vector 
is  changed  by  180  degrees  at  these  times.  Flux  is  a  product  of  the  current  and,  thus,  FOC 
experiences  increased  flux  when  the  current  is  increased  to  supply  the  demanded  torque. 

In  summary,  the  verification  of  V/f  control  through  the  comparison  of  measured 
and  simulated  results  allowed  FOC  and  DTC  to  be  simulated  with  confidence  that  they 
accurately  predict  a  physical  system  using  these  control  techniques.  These  three  control 
schemes  were  compared  to  one  another,  and  conclusions  are  discussed  in  the  next 
chapter. 
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V.  CONCLUSIONS 


A.  CONCLUSIONS 

Using  parameters  from  a  VSI-fed  induetion  maehine,  a  Simulink  model  of  this 
machine  was  designed,  simulated,  constructed  and  tested  using  standard  engineering 
principles  and  theory.  Thorough  research  and  simulation  phases  allowed  for  relatively 
short  development  and  construction  phases.  Experimental  results  demonstrated  that  this 
model  accurately  predicted  the  behavior  of  the  machine.  The  simulation  was  used  to 
evaluate  three  different  torque  control  schemes  (V/f,  FOC,  and  DTC). 

The  results  found  during  this  thesis  are  summarized  in  Table  12. 


Table  12.  Comparison  of  control  schemes,  from  [1]. 


Torque 

Flux 

Control  Type 

Con^x>l 

Control 

Response 

Advantages 

Disadvan  tages 

DC  Drive 


Direct 


Scalar  Frequency  None 

Control 


Flux  Vector  Control  Indirect 


Direct  Torque 
Control 


Direct 


Direct 


None 


Direct 


Direct 


High  High  accuracy 

Good  torque  response 
Simple 


Low  No  encoder 
Simple 


High  High  accuracy 

Good  torque  response 


High  No  encoder 

Moderate  accuracy 
Excellent  torque  response 


Motor  maintenance 
Motor  cost 

Encoder  required  for  high  accuracy 


Low  accuracy 
Poor  torque  response 


Encoder  always  required 


Encoder  required  for  high  accuracy 


The  evolution  from  DC  drives  to  various  forms  of  AC  drives  has  been  driven  by 
the  need  for  better  response  times,  higher  aeeuraey,  and  eontrol  of  torque  and  flux 
(direetly  or  indireetly).  Theory  and  prineiples  for  V/f  eontrol,  FOC,  and  DTC  were 
eovered  in  detail.  The  distinet  advantages  that  FOC  and  DTC  possess  over  V/f  eontrol 
allow  the  reader  to  eonelude  that  sealar  eontrol  is  not  a  preferred  eontrol  seheme.  FOC 
demonstrated  better  aeeuraey  mainly  beeause  of  the  ability  to  eontrol  the  eurrent,  but  this 
aeeuraey  eomes  at  the  neeessity  of  eurrent  sensors.  DTC  allows  for  direct  control  of  both 
the  torque  and  the  flux,  but  it  is  not  as  accurate  and  experiences  high  ripple  content  in 
both  the  flux  and  torque  components.  The  steady-state  current  of  DTC  was  significantly 


67 


lower  (less  electrieal  losses)  than  that  of  FOC.  If  the  maehine  is  going  to  be  used  mostly 
at  constant  speeds,  DTC  is  the  better  choice.  If  the  machine  changes  speeds  frequently, 
then  the  preferred  choice  is  FOC. 

With  a  reliable  model,  the  US  Navy  has  the  means  to  build  a  large  scale  prototype 
for  laboratory  purposes  and  future  research.  Future  research  on  this  topic  will  go  above 
and  beyond  current  publications  and  will  have  an  enabled  freedom  to  exploit  the 
capabilities  of  the  VSI-fed  induction  machine. 

B,  FUTURE  RESEARCH 

While  the  hardware  prototype  used  in  this  thesis  utilized  only  V/f  control,  a 
prototype  that  utilizes  the  other  two  control  schemes  can  be  built.  One  could  use  a  similar 
process  developed  here  to  verify  that  the  simulation  accurately  predicts  the  behavior  of  an 
induction  machine  using  FOC  and  DTC.  Improvements  leading  to  substantially  higher 
performance  of  both  control  schemes  should  be  investigated.  For  example,  current 
sensors  could  be  introduced  to  DTC  to  improve  accuracy  and  reduce  the  ripple 
component  of  torque  and  flux.  In  addition,  to  evaluate  whether  the  trends  developed  here 
are  consistent,  machines  of  different  sizes  should  be  studied.  Many  questions  remain,  but 
a  foundation  leading  to  those  answers  was  provided  by  this  thesis. 
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APPENDIX  A.  DATASHEETS 


Model  3231  Three-Phase  WooiidTIotor  [nducUon 
Motor 


Each  phase  of  the  stator  windings  of  this  motor  is  inde¬ 
pendently  terminated  and  identified  on  the  faceplate  to 
permit  operation  in  either  delta  or  star  (wye)  configura¬ 
tion.  The  rotor  windings  are  brought  out  to  the  faceplate 
via  external  slip  rings  ana  brushes.  This  machine  can  l3e 
used  as  a  wound-rotor  induction  naotor,  phase  shPer, 
single-phase  variable  coupling  transformer,  three-phase 
transformer,  selsyn  control,  frequency  converter  or 
asynchroftous  induction  generator.  The  speed  of  this 
machine  can  l)e  controlled  through  the  use  of  the  Three- 
Phase  Rheostat  (Model  8731). 
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SPECIFICATIONS 


Model  S211  OC  Motor/Gencfator 

120il208V-60Hz 

22W380V-S0HZ 

240/415  V- SO  Hz 

Power  Reqursfnent 

12Qi'208V 

220m)V 

240/415  V 

Ratng  Motor  Output  Power 

Generator  Ouput  Power 

Armature  Voltage 
Shint  Reid  Voltaje 
Ful  Load  Speec 

Ful  Load  Motor  Current 

Full  Load  Generator  Current 

175  W 

120  W 

now 

120W 

120V*DC 

220V-DC 

240  V- DC 

120  V -DC 

220V-DC 

240  V- DC 

1800  F/tnin 

1500  nmn 

1500  r^vi 

Z8A 

UA 

1.1  A 

1  A 

05  A 

0.5  A 

Physical  Characienstics  Dimensions  (H  x  W  x  0) 

NetWeigM 

308  x  201  x44C  mm  (12.1  X  II.Sx  17.3  h) 

14.1  kg  (Sib) 

Mo<iel  8221  Four-Pote  S<tuirTel-Cage  Induction  Motor 

12«208V-60Hz 

22O1G8O  V- 50  Hz 

240/415  V- SO  Hz 

Power  Requrement 

120/206  V 

220380 V 

240/415  V 

Ratng  Output  Power 

Stator  Voltage 

Full  Load  Speed 

FuN  Load  Cun«nt 

175  W 

120^18  V.  3i)hase 

220/380  V.3.phase 

240/415  V.  3-phase 

1670  nVnin 

1360  nhwr 

1305r/m»i 

I^A 

052  A 

0  46  A 

Ptiysxal  CharacJensncs  Dimensions  (H  x  W  x  0) 

Net  Weight 

308x291  X  440  mm  ( 12.1  X  11  5x17.3  in) 

135kg/297lbl 

Model  8231  Three-Phase  Wourwi-Rolor  Induction  Motor 

120008  V.  60  Hz 

220/380  V- SO  Hz 

240/415  V.  SO  Hz 

Power  Requrement 

i2o«oev 

220380 V 

2404157 

Ratng  Output  Power 

stator  Voltage 
Rotor  Voltage 

Full  Load  Speed 

Fiil  Load  Current 

175  W 

120/206  V.3phase 

220/380  V.Ophase 

240/415  V.Ofhase 

60/104  V.  3-phase 

1 10/100  V.3i>hase 

120/208  V.  3-phase 

1500  F/hiin 

1240  r^nti 

1315  r/nwi 

1.3A 

0.53  A 

0.48  A 

Physical  Charactehsdcs  Dimensiorts  (H  x  W  x  0) 

Net  Weight 

308  X  201  x440  mm  jtZI  x  11.5x  17.3  in) 

14  kg  (30.8  b) 
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POWER  SUPPLY 
MODEL  8821 


The  Power  Supply  provides  fixed  and  variable  AC 
and  DC  voltage  sources,  all  terminated  by  color-coded 
4  mm  safety  sockets.  Independent  circuit  breakers,  reset 
at  the  front  panel,  protect  the  input  to  and  output  from  the 
Power  Supply.  Indicator  lamps  monitor  the  presence  of 
input  voltage  in  each  phase.  When  a  phase  leg  of  the 
site's  power  service  is  out,  the  lamp  goes  off  to  reflect  this 
condition. 


A  voltmeter,  connected  through  a  selector  switch, 
monitors  the  variable  AC  and  DC  outputs  and  fixed  DC 
output  A  24  V  AC  output  provides  a  low-voltage  supply 
required  to  operate  other  EMS  equipment  such  as  meter¬ 
ing  modules  and  modules  used  in  the  Power  Eiectronics 
Training  System. 


SPECIFICATIONS 


Model  8321  -  Power  Supply 

1 20/208  V- 60  htz 

220/380  V- SO  H2 

240/41SV-50HZ 

Input  Line  Voltage 

Line  Current 

Service  knst^lation 

120/208  V 

220.'380V 

2401415  V 

15  A 

10  A 

20  A.  3'pha5e.  5  wres.  star  /Wye)-  conr>ected.  rtduding  neutral  and 
grouno 

Outputs  Faced  AC  3-Phase 

Variable  AC  3-Phase 

Variable  DC 

F'xed  DC 

Low  Power  AC 

I2aC08V-15A 

220380  7 -10  A 

240M15V-10A 

0-120/206  V-5  A 

0-22a380V-3A 

0-240M15V-3A 

0-120  V-8  A 

0-22GV-5A 

0-240  V-5A 

120V-2A 

220V-1  A 

240V-1  A 

24  V  -  3  A 

Wafl  Outlet  inchjded 

NEMAL21-20  1  NEMAL22-2D  |  CLIPSAL  56S0520 

Power  Cord 

3mnom 

Physical  Characteristics  Dimensions  (H  x  W  x  D) 

Net  Weight 

308  X  287  X  500  mm  (12.1  X  11-3  X  19.7  in) 

18  4  leg  (40.5  lb) 
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Laboratory  Setup 


MECHANICAL  LOADS 
MODELS  8911, 8913,  AND  8960 


Model  8911  Electrodynamometer 


The  Model  891 1  Electrodynamometer  is  equipped  with 
a  clear  plastic  faceplate  fitted  with  a  chrome  piano  hinge. 
It  can  be  coupled  to  a  0.2-kW  drive  motor  through  the 
use  of  the  Timing  Belt  (Model  8942).  The  faceplate  can 
be  lowered  to  provide  access  to  the  inside  for  coupling 
the  Electrodynamometer  to  the  drive  motor.  When 
closed,  the  'Septate  is  secured  by  two  quick-lock  fasten¬ 
ers. 

The  Electrodynamometer  is  made  of  a  squirrel-cage 
rotor  fitted  in  a  DC-excited  stator.  Mechanical  loading  is 
achieved  by  increasing  the  stator  field  current,  which 
generates  eddy  currents  in  the  driven  rotor.  The  stator 
is  trunnion-mounted  on  ball  bearings  and  rotates  against 
a  helicoidal  spring  to  oppose  braking  torque.  The  front- 
end  bell  is  provided  with  a  circular  scale  that  indicates  the 
torque  developed.  The  circular  scale  is  graduated  in 
either  Imperial  Units  (Ibfin)  or  Metric  Units  (N-m)  depend¬ 
ing  on  model  numbers.  The  Electrodynamometer  is 
electrically  powered  from  standard  fixed  AC  line  supply. 
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SPECIFICATIONS 


Model  8911  Electrodynamometer 

120/208  V-  60  Hz 

220/380  V  -  50  Hz  240/41 5  V  -  50  Hz 

Rating  Torque  Range  (891 1-1 ) 

Torque  Range  (8911) 

Speed 

Accuracy 

Input  Voltage 

Input  Current 

-0  3  to  +3  N  m 

-3  to  +27  Ibf  in 

250  to  3000  r/min 

2% 

120  V- AC 

220-240  V  -  AC 

2  A 

0.9  A 

Physical  Characteristics  Dimensions  (  H  x  W  x  D) 

Net  Weight 

308  x  291  X  490  mm  (12.1x11.5x19.3  in) 

17.4  kg  (38.3  lb) 

Model  8913  Prony  Brake 

Rating  Torque  Range  (8913-1 ) 

Torque  Range  (8913) 

Speed 

Accuracy 

0  to  3.4  N  m 

0  to  30  Ibf  in 

0  to  5000  r/min 

2% 

Physical  Characteristics  Dimensions  (  H  x  W  x  D) 

Net  Weight 

308  x  291  x440  mm  (12.1  X  11.5  X  17.3  in) 

8.8  kg  (19.4  lb) 

MECHANICAL  LOADS 
MODELS  8911,  8913,  AND  8960 


SPECIFICATIONS  (cont'd) 


Model  8960  Prime  Mover  /  Dynamometer 

120/208  V- 60  Hz 

220/380  V  -  50  Hz  240/41 5  V  -  50  Hz 

Prime  Mover  /  Dynamometer  Magnetic  Torque 

Rotation 

Display 

Torque  Units 

Friction  Torque 

Low  Power  Input 

Shaft  Encoder  Output 

Torque  Output  Signal 

Speed  Output  Signal 

Indicators 

0  to  3  N-m  (0  to  27  Ibf-in) 

cw  /  ccw 

Torque/Speed 

N-m  or  Ibf-in  (user  selectable) 

0.15  N-m  (1.3  Ibf-in)  @  1500  r/min 

24  V  -  AC 

A/B  -  360  pulses/r  -  TTL  Compatible 

0.3  N-nW  (2.655  Ibf-in/V) 

500  r/min  /V 

Overspeed,  Overheat,  Power  On 

Prime  Mover  Mechanical  Power 

Input  Voltage  Range 

Overvoltage  Detector 

Input  Current 

Speed  Range  (@  3.0  N-m  (26.6  Ibf-in)) 

750  W 

Oto  150  V-  DC 

NO 

YES 

6.5  A 

0  to  2500  r/min 

Dynamometer  Mechanical  Power 

Speed  Range  (@  3.0  N-m  (26.6  Ibf-in)) 

Speed  Range  (@  1.5  N  m  (13.3  Ibf  in)) 

Load  Control 

Load  Control  External  Input 

External  Load  Control  Input  Sensitivity 

250  W  continuous  -  750  W  maximum 

290  to  3000  r/min 

155  to  3000  r/min 

Manuel/Extemal 

Oto+10  V 

0.3  N  nW  (2.655  Ibf-in/V) 

Physical  Characteristics  Dimensions  (  H  x  W  x  D) 

Net  Weight 

308  x291  X  490  mm  (12.1  x  11.5  x  19.3  in) 

23.6  kg  (51.9  lb) 
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Manually 

Operated 

P/N 

Input 

Volts 

Phase 

Hertz 

* 

Output 

Volts 

Constant 

Current 

Rating 

I3PN2210B 

1  120  1 

1  1  PH  1 

|50/60  1 

|0-140  1 

|22A  1 

aPNtoooa 

3PNI2D0e 

3PNJC01BS«nM 

3PNisoae 

3PNJ401BS«m 

3PK221(ie 

3PN2S2t« 

mm  1  i.  rfihii..  . .  iiii>i-!-iWt-i  i..'  .ii 

mm  ..iii-i'.'iii-vti : .  ..ii 

W1  4«2-(1174fnn) 

624*115841111 

7  95*t2fl21rm) 

W7  3irtiaoanin» 

43r(1096mi 

4  31*  11096  im) 

K  1  06-(27Omi 

100*i254mn) 

165*  (41 8  mi 

B  4  3ri111  Itnru 

55(7(1397  mi 

7 17(181  Ontni 

111  37Sr(9S3ml 

4  75:(120  7  mi 

637(161  9imi 

IC  3SCrrW9mi 

587(1492ran)( 

F  02ri71iTm) 

Q37(11rm 

037  (8 1  rmu 

MN  029*174  tm)) 

031'(7lnni 

03r(95rnn) 
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APPENDIX  B.  MATLAB  M-FILES 


A,  INITIAL  CONDITIONS  FILE 

omega  =  2*pi*60;  %Modify  this  variable  to  change  reference  frame 

omega_b  =  2*pi*60; 

omega  in  =  2*pi*60; 

bypass=l ; 

if  bypass==l 

tstep=50e- 6 
tstop=4 ; 

else 

tstep=2e-6 
tstop= . 9 ; 

end; 

twopibyS  =  2*pi/3; 
poles  =  4; 

polesby2J  =  poles/2/ (2 . 04e-3*3) ; 

rs=l 6 ; 
rr=8 ; 

Xls=10.3569; 

Xm= 159.8850; 

Xlr=10.3569; 

Xm=51 . 6; 

rr=14;  %adjustment  to  rr 

rs=29;  %  increased  to  improve  model 

rsbyXls  =  rs/Xls; 
rrbyXlr  =  rr/Xlr; 

Xaq  =  1 / ( 1 /Xm+1 /Xls+1 /Xlr ) ; 

Xad  =  Xaq; 

XaqbyXls  =  Xaq/Xls; 

XaqbyXlr  =  Xaq/Xlr; 

XadbyXls  =  Xad/Xls; 

XadbyXlr  =  Xad/Xlr; 

V_phase  =  220*sqrt (2) /sqrt (3) ;  %Peak  value  used  in  Simulink  for  source 
modei 

XaqbyXm  =  Xaq/Xm; 

XadbyXm  =  Xad/Xm; 

%Steady  state  analysis 

omega_rotor  =  .79*omega  in; 

slip  =  (omega_in-omega_rotor) /omega_in; 

Z1  =  rs  +  j*omega  in/omega  b*Xls; 

Z2  =  j *omega_in/omega_b*Xm; 

Z3  =  rr/slip  +  j *omega_in/omega  b*Xlr; 
las  =  V_phase/sqrt (2) / (Z1  +  Z2*Z3/ (Z2  +  Z3) ) ; 
lar  =  -las* (Z2*Z3/ (Z2+Z3) ) /Z3; 

Te  =  3*poles/2 *Xm/omega  b*real ( j *con j ( las ) * lar ) ; 

Iqse  =  sqrt (2) *real (las) ; 

Idse  =  -sqrt (2) *imag (las) ; 

Iqre  =  sqrt (2 ) *real ( lar)  ; 
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Idre  =  -sqrt (2 ) *imag ( lar) ; 
psi_mqe  =  Xm* ( Iqse+Iqre) ; 
psi_mde  =  Xm* ( Idse+Idre) ; 
psi  qse  =  Iqse*Xls  +  psi_mqe; 
psi_dse  =  Idse*Xls  +  psi_mde; 
psi_qre  =  Iqre*Xlr  +  psi_mqe; 
psi_dre  =  Idre*Xlr  +  psi_mde; 

Te2  =  3/2*4/2* (psi_dse* Iqse  -  psi_qse*Idse) /omega_b; 

psi  mqic=psi_mqe; 

ps i_mdi c=ps i_mde ; 

psi_qsic=psi_qse; 

psi_dsic=psi_dse; 

psi_qric=psi_qre; 

psi  dric=psi_dre; 

omegar_ic  =  omega_rotor ; 

t_eventl= . 5 ; 
speedl=l *omega_b 
t_event2=2 ; 
speed2=0 . 8*omega_b 
t_event3=3 ; 
speed3=0 . 6*omega_b 

Vdc=300 ; 

Tsw=l/15000; 
gainl  =  1/2; 
gain2  =  1; 
oversample=2 ; 

Kp_i=25*4 ; 

Ki_i=500*6; 
swvO=  [0  0  0 ]  ; 
swvl= [1  0  0 ]  ; 
swv2=  [1  1  0 ]  ; 
swv3= [0  1  0 ]  ; 
swv4= [0  1  1 ]  ; 
swv5=  [0  0  1 ]  ; 
swv6=  [1  0  1 ]  ; 
swv7=  [1  1  1 ]  ; 


B.  MACHINE  PARAMETERS 

clc; 

clear  all; 

P  rated=746/4;  %  1/4  HP  machine 

wrm=1670/60*2*pi;  %rad/sec  mechanical 

T  rated=P  rated/wrm/ . 113  %in  in*lbf 
%measured  data 

11=  [.694  .770  .903  1.082  1.308]; 

12=  [.692  .761  .895  1.083  1.297]; 

Wl=  [87  121  160  203  249]  ; 

W2=  [44  13  23  60  93] ; 

speed=  [1789  1760  1726  1686  1640  0 . 1] /60*4*pi;  %radians  per  second 
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torque=  [0  3  6  9  12  18]; 
power=Wl+W2 ; 

%No  Load 

Vasnl=201/sqrt  (3) ; 

Iasnl= . 68 ; 
rs=11.91;  % 


XmXls  sq=  (Vasnl/Iasnl) ''2-rs^2; 

%  Eq. 

1 . 6 

in 

lab 

report 

P  core=  (138-3*Iasnl''2*rs)  ; 

%  Eq. 

1.4 

in 

lab 

report 

rc=3*Vasnl^2/P  core; 

%  Eq. 

1.5 

in 

lab 

report 

%Block;ed  rotor 
Vasbr=99/sqrt (3) ; 

Iasbr=2 .09; 

rr= (234/3/Iasbr^2-rs) ;  %  Eq.  1.9  in  lab  report 

Xls=sqrt  (  (Vasbr/Iasbr) '^2- (rs+rr) '^2 ) /2 ;  %  Eq.  1.10  in  lab  report 

XM=sqrt (XmXls_sq) -Xls  %  from  line  16  above 

Xlr=Xls; 

Vas=201/sqrt  (3) ; 

wb=2*pi*60;  %  base  omega,  page  152 

we=wb;  %this  means  we/wb=l 

rr_adj=8;  %adjustment  to  rr 
rs=16;  %  increased  to  improve  model 

wr= [0 : .5:0. 9999*wb] ; 
slip=(we-wr) ./we; 

Zl=rs  + j  *Xls ; 

Z2=rc*j*XM. / (rc+j*XM) ; 

Z3=rr_adj . / slip+ j  *Xlr; 

Zin=Zl  +  Z2  .  *Z3  .  / (Z2  +  Z3)  ; 

Ias=Vas . /Zin; 

Iar=-Ias .*Z2./(Z2+Z3); 

Ic=Ias  +  Iar . * ( j  *XM. / (rc+j  *XM)  )  ; 

Pcore=3*rc* Ic . *con j (Ic) ; 

Pscu=3*rs* las . *con j (las) ; 

Prcu=3*rr_adj  *Iar . *conj  (lar) ; 

Te=3*4/2*XM/wb*real ( j . *conj (las) .*Iar) ; 
power_out=Iar . *conj (lar) . *rr_adj . * ( 1-slip) . /slip; 

figure ( 1 ) ; 

plot (wr , Te/ 1 . 3558 17 9* 12 , ' linewidth ' , 2 ) ;  %torque  expressed  in  in*lbf 
hold  on; 

%  plot (wr, abs ( las ),' g linewidth 2 ) ;  %torque  expressed  in  in*lbf 
plot (speed, torque, ' ro ') ;  %torque  expressed  in  in*lbf 
hold  off; 

title (' Measured  and  Equivalent  circuit  torque  vs  speed'); 
xlabel ( ' Elec .  Speed,  rad.  per  sec.'); 

%  Stall  voltage,  section  8  in  the  lab 

speedb=[1722  1701  1665  1585  1713  1683  1629  1463] /60*4*pi; 
V_LLb=[199.4  176.8  157.1  136.5  187.3  166.8  146.9  126.3]; 
slipb= (we-speedb) ./we; 

Z3b=rr_adj . / slipb+ j  *Xlr; 

Zinb=Zl  +  Z2 . *Z3b. / (Z2  +  Z3b)  ; 
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Iasb=V_LLb. /sqrt (3) ./Zinb; 

Iarb=-Iasb. *Z2./ (Z2+Z3b) ; 

Teb=3*4/2*XM/wb*real (j.*conj (lasb) .*Iarb) /I. 3558179*12;  %torque 
expressed  in  in*lbf 
%  figure  (2 )  ; 

%  plot ( speedb, Teb, ' ro ' ) ; 

%  title (' Prediction  of  torque  for  the  stall  voltages  that  were 
measured ' )  ; 

%  xlabel ( ' Elec .  Speed,  rad.  per  sec.'); 


C.  SPEED  COMPARISON  FILE 

%Run  both  speed  simulations  first  to  create  the  data  for  these  plots 
close  all; 
clear  all; 


load  data_80_100b; 
la_80_100  =  add; 
lb_80_100  =  adc2; 
Vqs_80_100  =  adc3; 
RPM_80_100  =  adc4; 
load  data_100_80b; 
la_100_80  =  add; 
lb_100__80  =  adc2; 
Vqs_100_80  =  adc3; 
RPM_100_80  =  adc4; 
load  data_80_60b; 
la_80_60  =  add; 
lb_80_60  =  adc2; 
Vqs_80_60  =  adc3; 
RPM_80_60  =  adc4; 
time_shif tl= .44+. 0164; 
time_shif t2  =  l . 4+ . 0  6+ . 5 ; 
time  shift3=2 . 4+ . 055+ . 5; 


alpha=2*pi*500 ; 
lb_80_100_f il=zeros (1, 4095) ; 
lb_80_100_f il2  =  zeros (1, 4095)  ; 
lb_100_80_f il=zeros (1, 4095)  ; 
lb_100_80_f il2  =  zeros (1, 4095)  ; 
lb_80_60_f il=zeros (1, 4095)  ; 
lb_80_60_f il2  =  zeros (1, 4095)  ; 
for  ii=2 : 4095 

lb_80_100_f il (ii) =lb_80_100_f il (ii-1) +deltat*alpha* (lb_80_100 (ii 
1) -lb_80_100_fil (ii-1) )  ; 

lb_100_80_f il (ii) =lb_100_80_f il (ii-1) +deltat*alpha* (lb_100  80 (ii 
1) -lb_100_80_fil (ii-1) )  ; 

lb  80_60_fil (ii) =lb_80_60_f il (ii-1) +deltat*alpha* (lb_80_60 (ii-1) 
lb_80_60_fil (ii-1) )  ; 
end 

for  ii=2 : 4095 

lb_80_100_fil2 (ii) =lb_80_100_fil2 (ii- 
1) +deltat*alpha* ( lb_80_100_f il (ii-1) -lb_80_100_f il2 (ii-1) ) ; 
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Ib_100_80_fil2 (ii) =Ib_100_80_fil2 (ii- 
1) +deltat*alpha* (Ib_100_80_fil (ii-1) -Ib_100_80_f il2 (ii-1 
Ib_80_60_fil2 (ii) =Ib_80_60_fil2 (ii- 
1) +deltat*alpha* ( Ib_80_60_f il (ii-1) -Ib_80_60_fil2 (ii-1) ) 
end 

vbyf=0 ; 

DTC=1 

sim  ind_machine  3hp_speed  control_DTCver4 
simout_DTC=simout_vec; 
sw_state_DTC=sw_state; 
f luxout=f lux; 

vbyf=0 ; 

DTC=0 

sim  ind  machine_3hp_speed  control_DTCver4 
simout_FOC=simout_vec; 
sw_state_FOC=sw_state; 
f luxout_FOC=f lux; 

vbyf=l ; 

DTC=0 

sim  ind_machine_3hp_speed_control_DTCver4 
templ=length (simout_DTC) ; 
time_vec= [ 1 : tempi ] / tempi *tstop; 
sw_state_pwm=sw_state_l ; 

figure ( 1 ) ; 

plot ( time_vec, simout_vec ( : , 2 ) , ' g ' , ' linewidth ' , 3 ) ; 
hold  on; 

plot ( time+time_shif tl , RPM_80_100 , ' m ' , ' linewidth ' , 3 ) ; 
plot ( time+time_shif t2 , RPM_100_80 , ' m ' , ' linewidth ' , 3 ) ; 
plot ( time+time_shif t3 , RPM_80_60 , ' m ' , ' linewidth ' , 3 ) ; 
hold  off; 

title (' Measured  and  simulated  results  (V/f  control) '); 
legend ( ' V/f ' , ' Measured 'location', ' northeast ' ) ; 
xlabel('Time  (seconds)'); 
ylabel ( ' Speed  (rpms) ' ) ; 
grid  on; 

saveas  ( 1 ,  'figl.jpg',  'jpg'); 
figure  (2 ) ; 

plot ( time_vec, simout_vec ( : , 4 ) ,  '  g  '  ,  ' linewidth ' , 3 ) ; 
hold  on; 

plot ( time+time_shif tl , Ib_80_100_f il2 , ' m ' , ' linewidth ' , 3 ) ; 
%  plot (time+time_shiftl , Ib  80_100,'r'); 

plot ( time+time_shif t2 , Ib_100_80_f il2 , ' m ' , ' linewidth ' , 3 ) ; 
plot ( time+time_shif t3 , Ib_80_60_f il2 , ' m ' , ' linewidth ' , 3 ) ; 
hold  off; 

title (' Phase  B  current  filtered  at  300  Hz'); 

xlabel('Time  (seconds)'); 

ylabel ( 'Current  (amps) ' ) ; 

grid  on; 

if  vbyf  ==  1 


79 


legend ( ' V/f ' , ' Measured 'location', ' northeast ' ) ; 
else 

legend (' - ', ' - ', ' location ' , ' northeast ' ) ; 

end 

figure ( 3 ) ; 

plot  ( time_vec,  simout__vec  ( : ,  2 )  ,  '  g  '  ,  '  linewidth  '  ,  3 ) 
hold  on; 

plot ( time_vec, simout_FOC ( : , 2 ) , ' b ' , ' linewidth ' , 3 ) 
hold  off; 

title (' Speed  response:  V/f  vs.  FOC ' ) ; 
legend ( ' V/f ' , ' FOC ' , ' location ' , ' northeast ' ) ; 
xlabel('Time  (seconds)'); 
ylabel ( ' Speed  (rpms) ' ) ; 
grid  on; 

saveas  ( 1 ,  'figl3.jpg',  'jpg'); 
figure  ( 4 ) ; 

plot ( time_vec, simout_vec  ( : , 2 )  ,  '  g  '  ,  ' linewidth ' , 3 ) 
hold  on; 

plot ( time_vec, simout_DTC ( : , 2 ) , ' r ' , ' linewidth ' , 3 ) 
hold  off; 

title (' Speed  response:  V/f  vs.  DTC ' ) ; 
legend ( ' V/f ' , ' DTC ' , ' location ' , ' northeast ' ) ; 
grid  on; 

xlabel ( ' Time  (seconds)'); 
ylabel ( ' Speed  (rpms) ' ) ; 
saveas ( 1 , 'figl2.jpg', 'jpg'); 

figure ( 5 ) ; 

plot ( time_vec, simout_DTC ( : , 2 ) , ' r ' , ' linewidth ' , 3 ) 
hold  on; 

plot ( time_vec, simout_FOC ( : , 2 ) , ' b ' , ' linewidth ' , 3 ) 
hold  off; 

title (' Speed  response:  DTC  vs.  FOC'); 
legend ( ' DTC ' , ' FOC ' , ' location ' , ' northeast ' ) ; 
xlabel ('Time  (seconds)'); 
ylabel ( ' Speed  (rpms) ' ) ; 
grid  on; 

saveas  ( 1 ,  'figll.jpg',  'jpg'); 
figure  ( 6 ) ; 

plot (time_vec, simout_DTC  (  : , 4 ) ,  '  r  '  )  ; 
hold  on; 

plot (time_vec, simout_vec ( : , 4 ) , ' b ' ) ; 

title (' Phase  B  current  filtered  at  300  Hz'); 

legend ( ' DTC ' , ' FOC ' , ' location ' , ' southeast ' ) ; 

xlabel ('Time  (seconds)'); 

ylabel ( 'Current  (amps)  ' )  ; 

grid  on; 

figure  ( 7 )  ; 

plot (time_vec, simout_DTC ( : , 3) ,  ' r '  ,  ' linewidth ' , 3) 
hold  on; 


80 


plot ( time_vec, simout_FOC ( : , 3 ) , ' b ' , ' linewidth ' , 3 ) ; 
title (' Reference  Torque'); 

legend ( ' DTC ' , ' FOC ' , ' location ' , ' northeast ' ) ; 
xlabel ( ' Time  (seconds)'); 
ylabel (' Torque  (Newton  meters)'); 
grid  on; 

figure  (  8 )  ; 

plot ( time_vec, simout_DTC {:,!), ' r ' , ' linewidth ' , 3 ) ; 
hold  on; 

plot ( time_vec, simout_vec {:,!), ' b ' , ' linewidth ' , 3 ) ; 
%  axis([0  time_vec (tempi )  -2000  6000]); 

title (' Electrical  torque'); 

legend ( ' DTC ' , ' FOC ' , ' location ' , ' northeast ' ) ; 
xlabel ('Time  (seconds)'); 
ylabel (' Torque  (Newton  meters)'); 
grid  on; 


startl=29000; 
start2=30000 ; 


figure ( 9 ) ; 

plot ( f luxout (startl:start2, 1)  ,  f luxout (startl:start2,2)  ,  'r',  ' linewidth ' , 
3)  ; 

hold  on; 

plot ( f luxout_FOC (startl:start2, 1) , f luxout_FOC (startl:start2,2), 'b', 'lin 
ewidth ' , 3 )  ; 
hold  off; 

title (' Stator  flux  of  DTC  vs.  Stator  flux  of  FOC'); 
legend ( ' DTC ' ,  ' FOC '  ,  ' location '  ,  ' southeast ' ) ; 
xlabel (' Electrical  flux  (Volt  meters) '); 
ylabel (' Electrical  flux  (Volt  meters) '); 
grid  on; 

figure  ( 10 )  ; 

plot  ( time_vec,  sqrt  ( f  luxout  (  :  ,  1 )  .  2 +  f  luxout  (:,2)  .''2),  'r',  '  linewidth  '  ,  3 )  ; 
hold  on; 

plot ( time_vec, sqrt ( f luxout_FOC  (  : , 1 )  . ^2  +  f luxout_FOC (:,2)  .^2),  'b',  ' linewi 
dth ' , 3 )  ; 
hold  off; 

title (' Electrical  Flux'); 

legend ( ' DTC ' , ' FOC ' , ' location ' , ' northeast ' ) ; 
xlabel ('Time  (seconds)'); 
ylabel ('Flux  (Volt-meters)'); 
grid  on; 

sw_ct=0 ; 

for  ii=l:templ-l 

if  sw  state_DTC ( ii , 2 )  ==0  &  sw_state_DTC ( ii+1 , 2 )  ==  1 
sw_ct=sw_ct+l ; 

end 

end 
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sw_ctpwm=0 ; 
for  ii=l : templ-1 

if  sw_state_pwm ( ii , 1 )  ==0  &  sw_state_pwm ( ii+1 , 1 )  ==  1 
sw_ctpwm=sw_ctpwm+l ; 

end 

end 


D.  CHIPSCOPE  PLOTTER  FILE 

close  all; 
clear  all; 

datain=importdata ( ' im_datal00_80b.prn ' ) ; 
vecsize=length (datain . data) ; 
deltat=l/ (25e6/3600) ; 

%  datasize=round (1/60/deltat) ; 
datasize=vecsize-l; 
adclraw= zeros (l,datasize)  ; 
adc2raw= zeros (l,datasize) ; 
adc3raw=zeros (l,datasize)  ; 
adc4raw=zeros (l,datasize)  ; 
for  ii=l : datasize 

index=ii+vecsize-datasize- 1 ; 
if  datain . data ( index, 12+2 ) ==0 

add  raw  ( ii )  =datain  .data  ( index,  1+2)  +2*  da  tain,  data  ( index,  2+2)  +2'^2*datain. 
data ( index, 3+2 )+ . . . 

2^3*  da  tain,  data  ( index,  4+2)  +2'^4*  da  tain,  data  ( index,  5+2 )  +  .  .  . 

2^5* da tain. data ( index, 6+2) +2^6* da tain. data ( index, 7+2 )+ . . . 

2''7*  da  tain,  data  ( index,  8+2)  +2'^8*  da  tain,  data  ( index,  9+2 )+  .  .  . 

2^9* da tain. data ( index, 10+2) +2^10* da tain. data ( index,  11+2)  ; 

else 

add raw ( ii ) =datain .data ( index, 1+2) +2* da tain. data ( index, 2+2) +2^2* da tain, 
data ( index, 3+2 )+ . . . 

2^3*  da  tain,  data  ( index,  4+2)  +2'^4*  da  tain,  data  ( index,  5+2 )+  .  .  . 

2^5* da tain. data ( index, 6+2) +2^6* da tain. data ( index, 7+2 )+ . . . 

2^7*  da  tain,  data  ( index,  8+2)  +2'^8*  da  tain,  data  ( index,  9+2 )  +  .  .  . 

2''9*  da  tain,  data  ( index,  10+2)  +2^10*  da  tain,  data  ( index,  11+2)  -  2''11; 

end 

if  datain . data ( index, 12+14 ) ==0 

adc2raw ( ii ) =datain .data ( index, 1+14) +2* da tain. data ( index, 2+14) +2 ^2* data! 
n . data ( index, 3  +  14 )+ .  .  . 

2''3*  da  tain,  data  ( index,  4  +  14)  +2^4*  da  tain,  data  ( index,  5+14 )+  .  .  . 

2^5* da tain. data ( index, 6+14) +2^6* da tain. data ( index, 7+14 )+ . . . 

2''7*  da  tain,  data  ( index,  8  +  14)  +2^8*  da  tain,  data  ( index,  9+14 )+  .  .  . 

2^9* da tain. data ( index,  10  +  14) +2^10* da tain. data ( index,  11  +  14)  ; 

else 
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adc2raw  (ii)  =datain .  data  (index,  1  +  14)  +2*  da  tain,  data  (index,  2  +  14)  +2'^2*datai 
n . data ( index, 3  +  14 )+ .  .  . 

2^3* da tain. data ( index, 4+14) +2^4* da tain. data ( index, 5+14 )+ . . . 

2^5* da tain. data ( index, 6+14) +2^6* da tain. data ( index, 7+14 )+ . . . 

2^7* da tain. data ( index, 8  +  14) +2^8* da tain. data ( index,  9+14 )+ .  .  . 

2 ^ 9*datain . data  (index,  10  +  14)  +2'^10*datain.data  (index,  11  +  14)  - 

2^11; 

end 

if  datain . data ( index, 12+2 6 ) ==0 

adc3raw ( ii ) =datain .data ( index, 1+26) +2* da tain. data ( index, 2+26) +2^2*datai 
n . data ( index, 3+2 6 )+ . . . 

2^3* da tain. data ( index, 4+26) +2^4* da tain. data ( index, 5+26)+. . . 

2^5* da tain. data ( index, 6+26) +2^6* da tain. data ( index,  7+26)+.  .  . 

2''7*  da  tain,  data  ( index,  8+26)  +2^8*  da  tain,  data  ( index,  9+26)+.  .  . 

2^9* da tain. data ( index,  10+26) +2^10* da tain. data ( index,  11+26)  ; 

else 

adc3raw ( ii ) =datain .data ( index, 1+26) +2* da tain. data ( index, 2+26) +2^2*datai 
n . data ( index, 3+26)+. . . 

2^3* da tain. data ( index, 4+26) +2^4* da tain. data ( index, 5+26)+. . . 

2''5*  da  tain,  data  ( index,  6+26)  +2''6*  da  tain,  data  ( index,  7+26)+.  .  . 

2''7*  da  tain,  data  ( index,  8+26)  +2''8*  da  tain,  data  ( index,  9+26)+.  .  . 

2^9* da tain. data ( index,  10+26) +2^10* da tain. data ( index,  11+26)  - 

2''11; 

end 


adc4raw  (ii)  =da  tain,  data  (index,  1  +  38)  +2*  da  tain,  data  (index,  2  +  38)  +2'^2*datai 
n . data ( index, 3  +  38 )+ . .  . 

2^3* da tain. data ( index, 4+38) +2^4* da tain. data ( index, 5+38 )+ . . . 

2^5* da tain. data ( index, 6+38) +2^6* da tain. data ( index, 7+38 )+ . . . 

2^7* da tain. data ( index, 8  +  38) +2^8* da tain. data ( index,  9+38 )+  .  .  . 

2 9*datain . data  (index,  10  +  38)  +2'^10*datain.data  (index,  11  +  38)  ; 


end 

adcl=adclraw/2  6 ; 
adc2=adc2raw/2^  6 ; 
adc3=adc3raw/2’^ll*251/sqrt  (3)  ; 
adc4=adc4raw; 

time= [0 :datasize-l] *deltat; 

set ( 0 , ' Def aultAxesFontName ' , ' times ' ) 
set(0,  ' DefaultAxesFontSize '  ,  14 ) 
set  ( 0 ,  ' Def ault Text FontName ' ,  ' times ' ) 
figure ( 1 ) ; 

plot  (time,  add ,  '  b  '  ,  '  linewidth  '  ,  2 )  ; 
hold  on; 

plot (time, adc2 , ' g ' , ' linewidth ' , 2 ) ; 
hold  off; 
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xlabel ( ' Time ' ) 
ylabel('DC  bus  Magnitude') 
grid  on 

figure  (2 )  ; 

plot (time, adc3 , ' r ' , ' linewidth ' , 2 ) ; 
hold  on; 

plot (time, adc4 , ' m ' , ' linewidth ' , 2 ) ; 

hold  off; 

xlabel ( ' Time ' ) 

ylabel ( ' Magnitude ' ) 

grid  on 

%axis([.0055  .006  -1  1]); 
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APPENDIX  C.  SIMULINK  MODEL 


A,  INDUCTION  MACHINE  MODEL 


Induction  Machine  model  A  Julian.  Mar  2013 
Equations:  4.14-1  to  4.14-19  in  Krause  text 
1/4  HP  parameters  from  lab  machine 
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B 


CONTROLLER 


1,  Inverter 
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Transformation 


Reference  1 


^  theta  Sectors  j- 
Sector  selectar 


DotProductI 


DotProduct2 


c.  Sector  Selector 


Q> 

thelB 


pi/3 


2*pi/3-f 


-pi/3 


■2-pi/3  J  ^ 


- ^ 


-KZ) 


d.  Subsystem 
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Subsystem  1 


L 

TOM 


Zone  1-5  Test 


g.  Subsystem  2 


Inssgef  jnserlEO  here  is  ttie  (sampler  frequencyViPWM  frequency! 


h.  Triggered  Subsystem 


Scope2 
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1. 


Voltage 


Add1 


Inverse  Transformation 


Inverse  K_  Transformation 


c 


INDUCTION  MACHINE 


1.  DTC 
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First  Row 


h 


Second  Row 


Indec 

Vect<j1 
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Third  Row 


n 


d. 


Fourth  Row 
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f.  Sixth  Row 


3. 


(Pdr 


4.  (Pds 


Vmd 
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XaqbyXIr 


9. 


Rotor  Flux 


fqsx 


10,  Stator  Flux 


fqsx 


11,  Torque 


Changes  electrical  torque  from  perunit  to  absolute  values 
Eq  4  6-5 
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D. 


SPEED 


electrical  radians  per  second 


mechanical  revolutions  per  minute 


Gaini 


KjD 

Outi 


E.  TORQUE  EQUATION 


F.  INVERSE  TRANSFORMATION 


G.  LOAD  TORQUE 


In1 


u2 

J 

- ► 

Math 

Function 


-K- 
Gain 


-KjD 

Outi 
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